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SURFACE  CRACKS  EMANATING  FROM  FASTENER  HOLES 


The  finite  element-alternating  method  is  used  to  determine  stress 
intensity  factors  along  the  periphery  of  a part-elliptical  crack 
emanating  from  a fastener  hole  in  a finite-thickness  plate.  A sequence 
of  iterations  is  performed  between  an  analytic  solution  for  an  ellipti- 
cal crack  embedded  in  an  infinite  solid  and  a finite  element  solufon 
for  the  uncracked  plate  with  a fastener  hole.  The  technique  determines 
the  stress  field  near  the  crack,  the  stress  intensity  factor  and  the 
crack  opening  displacements. 

Mode-one  stress  intensity  factors  around  the  crack  front  are 
presented  for  three  classes  of  crack  location  relative  to  the  hole  and 
numerous  crcivk  shapes  and  sizes.  Calculations  are  performed  for  cracks 
emanating  from  both  loaded  and  unloaded  fastener  holes.  Crack  opening 
displacements  for  all  cases  are  presented. 

The  results  of  this  study  are  compared  wlth^static  fracture  tests 
in  polymethylmethacrylate  and  with  experiments  and  estimates  of  other 
authors. 
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CHAPTER  1 


INTRODUCTION 

1 .1  Motivation  for  Study  of  Cracks  Emanating  from  Fastener  Holes 

Cracks  in  aerospace  structures  have  caused  catstrophic  failures 
resulting  in  loss  of  life,  destruction  of  millions  of  dollars  in  equip- 
ment and  severely  reduced  operational  capabilities.  While  the  elimina- 
tion of  possibilities  for  crack  initiation  and  growth  is  a worthwhile 
endeavor,  the  fact  remains  that  flaws  due  to  material  defects,  manu- 
facturing methods,  and  in-service  conditions  will  always  be  present  to 
some  extent.  This  simple  admission  opens  the  door  to  design  against 
failure  in  the  presence  of  flaws  existing  at  the  onset  of  operational 
use  and  during  periods  of  fatigue  growth.  The  fact  that  a crack  exists 
is,  by  Itself,  no  longer  a criteria  for  the  scrapping  of  an  aerospace 
system  nor  is  it  sufficient  reason  for  extensive  modifications. 

The  United  States  Air  Force  is  currently  developing  a design 
criteria  (1)*  to  reduce  the  possibility  of  catastrophic  failures  caused 
by  crack-like  defects.  This  criteria  will  be  established  on  the  basis 
of  three  separate  yet  interdependent  studies: 

1.  At  the  beginning  of  a particular  operational  period,  determine 
the  most  serious  initial  crack  which  could  exist  in  a struc- 
ture and  predict  the  most  severe  operational  conditions  that 
could  occur. 

2.  Determine  the  crack  sizes  which  would  cause  catastrophic 
failure. 

♦Numbers  in  parentheses  refer  to  references  listed  at  the  end  of  the 
text. 


1 


2 


3.  Using  predicted  system  operational  data,  calculate  the  time  or 
number  of  cycles  for  existing  cracks  to  grow  to  critical  size. 
The  first  and  last  studies  listed  are  largely  measurement  oriented  and 
thus  perhaps  more  easily  completed  than  the  frequently  complicated 
theoretical  problem  posed  by  the  second  study.  It  is  to  the  partial 
completion  of  the  second  study  that  this  work  is  dedicated. 

A comprehensive  review  of  aerospace  structural  failures  completed 
by  the  United  States  Air  Force  in  1971  (2)  showed  the  origin  of 
failures  due  to  cracks,  in  order  of  increasing  frequency  of  occurrence, 
to  be: 

1 . Surface  cracks 

2.  Corner  cracks 

3.  Cracks  emanating  from  fastener  holes 

A full  one-third  of  all  failures  studied  were  due  to  cracks  emanating 
from  fastener  holes.  While  progress  in  the  solution  of  surface  and 
corner  crack  problems  has  been  promising,  only  two  direct  three- 
dimensional  solutions  (3,4)  are  known  to  date  for  the  three-dimensional 
problem  of  cracks  emanating  from  fastener  holes.  This  work  presents 
stress  intensity  factor  solutions  for  three  crack  orientations  at  open 
and  loaded  fastener  holes  in  a finite-thickness  plate. 

1.2  Previous  Work  Leading  to  Stress  Intensity  Factors  for  Cracks 

Emanating  from  Open  Fastener  Holes 

The  first  major  advance  in  the  study  of  cracks  emanating  from  open 
holes  in  plates  was  due  to  Bowie  (5)  in  1956.  He  solved  the  two  dimen- 
sional problem  of  single  and  double  through-cracks  at  an  open  hole 
(Figures  1.1  and  1.2)  using  complex  variable  methods.  Although  not  an 
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exact  solution,  Bowie's  results  have  served  as  a basis  in  the  framing 
of  other  studies.  Perhaps  more  important  is  the  almost  universal  use 
of  his  solution  to  establish  the  accuracy  of  other  methods  designed  to 
provide  answers  for  the  more  complicated  problems  being  attempted  today. 

Tweed  and  Rooke  (6)  considered  the  single  through-crack  problem 
at  a later  date.  They  used  a Mel  1 in  transform  technique  to  predict 
stress  intensity  factor  relationships.  Their  solution  is  claimed  to 
have  greater  accuracy,  especially  for  smaller  crack  lengths,  than  the 
Bowie  solution. 

Since  the  publication  of  Bowie's  work,  theoretical  approaches  to 
the  more  complicated  three-dimensional  problem  of  a surface  flaw 
emanating  from  a fastener  hole  in  a plate  have  been  limited  largely 
to  estimates  based  upon  extensions  of  related  problem  solutions.  In 
general,  these  estimates  suffer  from  an  incomplete  application  of  the 
many  factors  influencing  the  problem  and  from  oversimplifications. 
Nevertheless,  it  is  worthwhile  to  review  some  of  these  estimates  to 
impress  upon  the  reader  the  various  directions  taken  in  the  study  of 
this  problem. 

Kobayashi  (7)  estimated  stress  intensity  factors  for  a semi- 
elliptical embedded  crack  adjacent  to  an  open  fastener  hole  in  a very 
thick  plate  (Figure  2.1c).  A linear  approximation  to  the  plane  strain 
stress  distribution  near  the  hole  was  applied  to  a circular  crack  in 
an  infinite  solid  without  a hole.  The  stress  intensity  factor  for  an 
elliptical  crack  was  formed  through  an  estimated  shape  correction  to 
the  results  for  the  circular  crack.  Surface  effects  and  the  through- 
thickness  stress  variation  seen  in  the  three-dimensional  problem  are 
neglected. 
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Hsu  and  Liu  (8)  estimated  stress  intensity  factors  for  a corner 
crack  at  a hole  in  a plate  (Figure  2.1a)  by  using  as  a basis  the  solu- 
tion for  an  elliptical  crack  in  an  infinite  solid.  They  applied  a 
series  of  estimated  correction  factors  to  account  for  the  presence  of 
the  hole  and  plate  surfaces.  Bowie's  two-dimensional  through-crack 
solution  was  assumed  to  provide  a valid  hole  correction  while  an  inter- 
polation between  three  special  crack  shapes  produced  a front  surface 
correction.  No  back  surface  correction  was  introduced. 

Liu  (9)  considered  a quarter-circular  crack  at  a hole  in  a plate. 

He  applied  Smith's  (10)  solution  to  approximately  account  for  the  hole 
surface  and  front  surface,  Kobayashi's  (11)  solution  as  an  approximate 
back  surface  correction,  and  Bowie's  (5)  two-dimensional  solution  to 
approximate  the  three-dimensional  hole  effect.  The  resulting  stress 
intensity  estimate  was  limited  to  a point  on  the  crack  periphery  midway 
between  the  front  and  hole  surfaces. 

Perhaps  one  of  the  more  complete  estimates  is  that  of  Shah  (12). 

He  began  by  applying  the  plane  strain  stress  distribution  near  a hole 
in  a plcite  under  uniaxial  loading  to  a circular  crack  and  used  a Green's 
function  to  find  stress  intensity  factors.  He  reasoned  that  these  re- 
sults also  applied  to  two  embedded  semi-circular  cracks  emanating  from 
a fastener  hole.  Estimated  crack  shape  correction  factors  were  applied 
to  extend  the  estimate  to  elliptical  shapes  and  Shah  and  Kobayashi's 
(13)  results  were  applied  as  approximate  back  surface  correction  fac- 
tors. A constant  front  surface  correction,  independent  of  location  on 
the  crack  border,  was  assumed.  A final  factor  was  applied  to  relate 
single  crack  results  to  those  for  double  cracks.  Shah  also  used  his 
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I procedure  to  estimate  stress  intensity  factors  for  cases  of  cracks 

I emanating  from  a neat-filled  loaded  fastener  hole. 

Grandt  (14)  recently  applied  a linear  superposition  technique 
I suggested  by  Rice  (15)  to  the  Bowie  problem.  He  extended  Bowie's  re- 

sults to  cases  of  through-cracks  adjacent  to  cold  worked  holes,  open 
I holes  with  varying  remote  loads,  and  interference-fit  loaded  fasteners. 

I A necessary  prerequisite  for  this  method  1s  a knowledge  of  stress  in- 

tensity factors  and  crack  opening  displacements  for  one  particular 
I crack  shape  and  orientation.  The  current  lack  of  such  information  has 

limited  extension  of  Grandt's  results  beyond  the  two-dimensional  through- 
I crack  problem. 

a The  very  nature  of  the  estimates  of  Kobayashi  (7),  Hsu  and  Liu 

(8),  Liu  (9),  Shah  (12)  and  others  has  caused  many  investigators  to  in- 
I tensify  their  efforts  in  seeking  results  through  experimental  means. 

Hall  and  Finger  (16)  completed  a series  of  static  fracture  tests 
I on  specimens  fashioned  from  2219-T87  aluminum  and  5A1-2.5  Sn  (ELI) 

I titanium  alloys.  They  suggested  two  empirical  failure  criteria  for  one 

location  on  the  crack  front  and  restricted  to  the  crack  geometries 
I tested. 

Hall,  Shah,  and  Engstrom  (17)  performed  similar  experiments  in 
I 4340  steel  plates  along  with  a wide  range  of  fatigue  and  spectrum 

m loading  tests. 

Static  fracture  tests  of  engineering  metals,  while  useful  when 
f specific  data  on  a particular  material  is  needed,  are  subject  to  diffi- 

culties when  compared  to  theories  based  upon  principles  of  linear 
^ elastic  fracture  mechanics.  If  specimen  sizing  is  too  small,  gross 

~ plastic  deformations  can  produce  results  which  are  quite  different  from 
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elastic  predictions.  Additionally,  the  point  of  initial  failure  on 
the  crack  border  is  difficult  to  locate  in  metals.  This  has  led  many 
investigators  to  the  use  of  transparent,  brittle  polymers  such  as 
polycarbonate,  epoxy,  and  polymethylmethacrylate  (PMMA)  as  better 
materials  for  modeling  the  fracture  process. 

Jolles,  McGowan,  and  Smith  (18)  experimentally  studied  the  single 
crack  adjacent  to  an  open  fastener  hole  using  established  stress 
freezing  photoelasticity  techniques.  They  extended  their  earlier  work 
(19)  to  more  crack  shapes  and  added  the  calculation  of  stress  intensity 
factors  at  two  intermediate  locations  on  the  crack  front  to  those  at 
the  hole  and  surface  intersections. 

Snow  (20)  followed  the  lead  of  Grant  and  Hinnerichs  (21)  with  a 
series  of  fatigue  experiments  in  PMMA.  Single  cracks  adjacent  to  open 
fastener  holes  were  introduced  and  fatigued  while  crack  growth  data  was 
recorded  photographically.  Suitable  baseline  tests  produced  additional 
information  allowing  for  solution  of  the  Paris  fatigue  relationship  for 
stress  intensity  factors.  Results  are  published  only  for  the  inter- 
sections of  the  crack  front  with  the  hole  and  with  the  front  surface. 

While  the  techniques  of  experimental  fracture  mechanics  are  surely 
becoming  more  sophisticated  and  capable  of  producing  some  credible 
results,  theoretical  solutions  to  the  three-dimensional  problem  are 
still  few  in  number. 

Browning  (3)  formed  a solution  based  upon  the  Schwartz-Neumann 
alternating  method  to  obtain  results  for  a single,  embedded,  semi- 
circular crack  emanating  from  an  open  fastener  hole  in  a finite  thick- 
ness plate  (Figure  2.1c).  Two  solutions  were  used  in  the  alternating 
procedure: 
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1.  Stress  in  a finite  body  due  to  surface  tractions  using  a 
three-dimensional  finite  element  method  (Figure  1.3). 

2.  Stress  in  an  infinite  body  due  to  arbitrary  normal  loading 
applied  to  an  embedded  circular  crack. 

Iteration  between  these  two  solutions  removes  tractions  on  the  finite 
body  free  surfaces  and  the  crack  surface.  The  accumulated  residuals 
from  freeing  of  the  crack  surface  are  used  to  compute  stress  intensity 
factors  from  the  circular  crack  solution.  Browning's  results  showed 
trends  similar  to  the  estimates  of  Kobayashi  (7).  Unfortunately,  his 
solution  took  a large  amount  of  computer  time  and  suffered  some  diffi- 
culties in  the  interpolation  of  stresses  from  finite  element  nodal 
points  to  Gaussian  quadrature  points.  This  necessitated  a specific 
finite  element  mesh  arrangement  for  each  crack  shape  considered. 
Nevertheless,  Browning's  pioneering  work  provided  the  basis  for  further 
studies. 

Ganong  (4)  extended  the  work  of  Browning  by  replacing  solution  2 
indicated  above  with  the  solution  of  Shah  and  Kobayashi  (22)  for 
stresses  in  an  infinite  medium  due  to  arbitrary  normal  loading  on  an 
embedded  elliptical  crack  (Figure  1.4).  He  also  decreased  the  compu- 
tational time  required  by  using  Cholesky  decomposition  of  the  finite 
element  stiffness  matrix  and  direct  calculation  of  stresses  at  the 
Gaussian  quadrature  points  on  the  finite  element  models.  Ganong  pre- 
sented results  for  a number  of  single  quarter-elliptical  cracks  emana- 
ting from  open  fastener  holes.  The  versatility  of  the  finite  element- 
alternating  method  has  several  facets: 

1.  Once  the  finite  element  stiffness  matrix  is  decomposed,  stress 
intensity  factor  results  for  a wide  variety  of  crack  shapes 
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can  be  found  at  the  expense  of  only  250  seconds  of  COC  6400 
central  processor  time  per  crack  problem. 

2.  The  finite  element  portion  of  the  alternating  method  controls 
the  geometry  of  the  finite  body  containing  the  crack  and  is 
capable  of  modeling  the  complex  geometries  of  actual  machine 
parts. 

3.  The  procedure  gives  stress  intensity  factors  directly  and 
does  not  rely  on  extrapolating  stresses  or  fitting  crack 
opening  displacements  as  do  some  of  the  direct  finite  element 
schemes. 

1 .3  Previous  Work  on  Stress  Fields  in  Plates  with  Loaded  Fastener  Holes 

The  first  step  of  the  finite  element-alternating  method  applies  a 
known  pressure  distribution  to  the  crack  plane.  In  the  case  of  an  open 
hole,  the  three-dimensional  analytic  solution  of  Sternberg  and  Sadowsky 
(23)  is  used.  Stress  intensity  factor  results  for  certain  cases  of 
load  transfer  to  a filled  fastener  hole  can  be  found  in  the  same  manner 
if  the  crack  plane  stress  distribution  is  known.  The  necessity  for 
knowledge  of  crack  plane  stress  distributions  is  not  unique  to  the 
finite  element-alternating  method.  Most  of  the  estimates  discussed 
earlier  have  the  same  requirement.  While  both  two-dimensional  and 
three-dimensional  stress  distributions  are  available  for  the  open  hole 
case,  the  complexity  of  the  loaded  fastener  hole  problem  has  limited 
available  results  to  two-dimensions  and  less  generality.  Figure  1.5 
shows  the  pin-loaded  hole  case  considered  in  this  work  where  the  ratio 
Q/P  ranges  from  0.0  for  an  open  hole  to  1.0  for  100%  load  transfer  to 
the  fastener. 
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Perhaps  the  earliest  attempt  to  define  the  stress  distribution 
near  a pin-loaded  hole  in  a plate  was  a series  of  photoelastic  studies 
by  Coker  (24)  during  the  years  1913  to  1919.  He  was  the  first  to  notice 
very  large  tensile  stresses  on  the  plane  of  minimum  cross-section  in 
a plate  of  limited  width. 

Bickley  (25),  at  the  urging  of  Coker,  used  a classical  elasticity 
approach  in  two  dimensions  to  analytically  find  stresses  around  a 
circular  hole  in  an  infinite  elastic  solid  due  to  surface  tractions  on 
the  hole  boundary.  Results  for  several  realistic  hole  surface  tractions 
were  calculated. 

In  1935,  Knight  (26)  made  a significant  contribution  with  a 
theoretical  solution  which,  even  today,  has  not  been  greatly  improved 
upon.  He  recognized  that  hole  surface  tractions  were  important  and 
assumed  their  form  from  a study  by  Howland  and  Stevenson  (27).  Using 
general  elasticity  theory,  Knight  then  accounted  for  a finite  plate 
width  of  twice  the  hole  diameter.  Theocaris  (28)  extended  Knight's 
work  to  three  additional  hole  diameter  to  plate  width  ratios  in  1956. 

Twelve  years  later,  Gregory  (29)  solved  the  same  problem 
(Figure  1.5)  with  the  theory  of  infinite  systems  of  linear  equations. 

For  tractions  on  the  hole  boundary,  he  assumed  that 


-3/4  (^)  (1-COS2B) 

d 


rr 


for  0 B £ It 
for  -Tf  < B < 0 


(1.3.1) 


and 


Org  =0  for  all  B 


(1.3.2) 


where  a is  the  hole  radius,  h the  plate  half-width,  T the  remote 


•J*. . 


10 


uniaxial  tensile  stress,  and  6 the  angle  measured  from  the  hole  edge  at 
the  minimum  cross-section.  This  boundary  condition  simulates  a smooth, 
rigid,  neat-fit  pin.  Gregory  presented  only  stress  concentration 
factor  results  at  6=0  for  a large  range  of  hole  diameter  to  plate  width 
ratios,  a/h.  These  results  compared  well  with  experiments,  particularly 
at  the  lower  ratios  of  a/h. 

The  work  of  Bickley  (25),  Knight  (26),  Theocaris  (28),  and  Gregory 
(29)  represent  the  extent  of  present  analytic  solutions  for  the  pin- 
loaded  hole  problem.  Many  important  effects  remain  unaccounted  for  in 
their  work  and  the  increased  problem  complexities  cause  more  recent  in- 
vestigations to  be  directed  toward  experimentation  and  approximation 
with  finite  element  methods. 

Frocht  and  Hill  (30)  used  photoelasticity  with  bakelite  specimens 
and  strain  gage  mountings  on  aluminum  plates  to  experimentally  de- 
termine stress  concentration  factors  for  plates  with  pin-loaded  holes. 

In  addition  to  a wide  range  of  plate  geometries,  they  tested  effects  of 
differing  pin  clearances  and  lubrications. 

Jessop,  Snell,  and  Hollister  (31)  also  applied  photoelastic  tech- 
niques to  the  problem,  considering  the  interaction  of  various  pin  load 
levels  with  three  different  pin  clearances.  Again,  results  were 
limited  to  the  pin-loaded  hole  boundary. 

Cox  and  Brown  (32)  reviewed  and  coordinated  available  experimental 
data  in  1963  in  an  attempt  to  sort  out  and  explain  the  many  results. 

In  addition,  they  clarified  the  theoretical  background  as  it  applied 
to  the  various  findings. 

These  studies,  while  enlightening,  consider  only  one  or  two  special 
effects  of  pin-plate  interaction  and  generally  present  rather  limited 
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results.  The  increasing  flexibility  of  the  finite  element  method  lends 
itself  to  studies  of  the  many  effects  acting  simultaneously. 

Harris,  Ojalvo,  and  Hooson  (33)  solved  a wide  variety  of  pin-load 
problems  using  finite  element  methods.  They  considered  both  linear 
and  non-linear  stress  and  deformation  of  mechanically  fastened  joints 
and  included  effects  of  pin  clearance  and  interference,  load  levels, 
plate  and  pin  geometry,  residual  stresses,  and  fastener  bending. 

Mohaghegh  (34)  prepared  a report  including  22  combinations  of 
fastener  and  plate  loading,  again  with  finite  element  methods.  Two- 
dimensional  results  are  presented  in  the  form  of  stresses  around  the 
periphery  of  a fastener  hole  containing  a rigid,  smooth  pin. 

Brombolich  (35)  considered  the  effects  of  plasticity,  load 
sequence,  and  fastener  interference  and  contact.  His  results  again 
were  obtained  using  finite  element  methods  and  are  restricted  to  un- 
loaded fasteners  in  plates  under  remote  uniaxial  tension. 

Due  to  the  many  effects  present  in  the  case  of  a pin-loaded  hole 
in  a plate,  choice  of  a suitable  crack  plane  stress  distribution  for 
adaptation  to  a stress  intensity  factor  solution  is  not  as  straight- 
forward as  it  might  at  first  seem.  In  all  likelihood,  an  appropriate 
distribution  for  a specific  case  understandably  will  not  be  found  in 
the  literature  and  it  becomes  necessary  to  generate  needed  results 
independently. 

1 .4  Objectives 

The  objectives  of  this  work  were  to: 

1.  Extend  Ganong's  (4)  analysis  to  more  crack  shapes  and  to 

two  additional  classes  of  crack  location  relative  to  the  hole. 
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2.  Apply  the  finite  element -alternating  method  to  cases  of  cracks 
emanating  from  pin-loaded  fastener  holes. 

3.  Design  and  perform  a series  of  static  fracture  tests  to  verify 
the  theoretical  results. 

4.  Calculate  mode-one  crack  opening  displacements  for  cracks 
emanating  from  fastener  holes. 

5.  Assess  accuracy  by  comparing  the  theory  with  experiments  of 
other  investigators. 

In  the  process  of  fulfilling  these  objectives,  the  finite  element- 
alternating  method  will  be  shown  to  accurately  and  efficiently  calculate 
stress  intensity  factors  for  a variety  of  cracks  emanating  from  both 
open  and  loaded  fastener  holes  in  finite-thickness  plates. 


CHAPTER  2 


THEORETICAL  APPROACH 

2.1  General  Problem  Statement 

Consider  part-elliptical  cracks  emanating  from  a hole  in  a finite 
thickness  plate  (Figure  2.1).  The  crack  center  is  located  at  the  hole 
edge  and  the  crack  lies  in  a plane  perpendicular  to  the  axis  of  loading. 
The  finite  plate  geometry  and  loading  conditions  are  such  that  the 
plane  of  the  crack  is  a plane  of  symmetry  for  open  hole  problems. 
Mode-one  stress  intensity  factors  are  to  be  calculated  around  the  crack 
border  for  cases  of  remote  uniaxial  loading  on  a plate  with  an  open 
fastener  hole,  and  partial  and  total  transfer  of  load  to  a filled 
fastener  hole  (Figure  1.5). 

2.2  General  Solution  Method 

The  finite  element-alternating  method  used  in  this  work  separates 
the  problem  into  two  problems  with  stress  intensity  factor  results 
obtained  through  an  iterative  procedure.  Figure  2.2  schematically 
presents  two  iterations  of  the  procedure  where  the  fastener  hole  has 
been  omitted  in  the  interest  of  clarity.  The  top  row  of  figures 
represents  Solution  1 while  the  bottom  row  shows  Solution  2.  In  the 
description  to  follow,  letters  in  parentheses  refer  to  steps  in  the 
iteration  process  as  shown  on  the  figure. 

The  first  solution,  hereafter  referred  to  as  Solution  1,  is  that 
of  stresses  in  an  unflawed  finite  body  due  to  the  applied  remote 
loading.  Solution  1 produces  a normal  stress  distribution,  p^^,  on  the 
plane  where  the  crack  is  to  be  located  due  to  remote  loads  (Step  a). 

The  second  solution,  hereafter  referred  to  as  Solution  2,  is  that  of 
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stresses  in  an  infinite  body  due  to  a pressure,  -p^,  applied  to  the 
surface  of  an  embedded  elliptical  crack.  The  solution  to  this  problem 
results  in  stresses  at  locations  in  the  infinite  body  corresponding  to 
points  on  the  surface  of  the  finite  body  and  frees  the  crack  of  residual 
stress  (Step  b).  The  stresses  calculated  by  Solution  2 at  surface  lo- 
cations on  the  finite  body  are  reversed  to  free  its  surfaces  and  Solu- 
tion 1 again  produces  a normal  pressure  distribution  at  the  location  of 
the  crack  in  the  unflawed  finite  body  (Step  c).  Iteration  between  the 
two  solutions  is  continued  until  stresses  on  the  finite  body  surfaces 
predicted  by  Solution  2 and  stresses  at  the  crack  location  from  Solu- 
tion 1 are  small.  The  superposition  of  results  from  all  iterations 
then  gives  a close  approximation  to  the  desired  solution  for  remote 
uniaxial  loading  on  a finite  body  with  a crack  where  the  crack  surface 
is  stress  free  and  the  finite  body  surfaces  have  the  prescribed  stress 
boundary  conditions.  The  total  of  normal  stresses  applied  to  the  crack 
from  all  iterations  to  free  its  surfaces  produces  mode-one  stress 
intensity  factors,  Kj,  which  are  computed  directly  from  Solution  2. 

It  should  be  mentioned  that  in  actual  practice.  Solution  1 is  not 
used  to  define  the  initial  distribution  of  pressure  on  the  crack  plane 
for  remote  uniaxial  loading  (Step  a).  Rather,  the  three-dimensional 
analytic  solution  of  Sternberg  and  Sadowsky  (23)  for  the  same  loading 
case  provides  the  necessary  stress  distribution. 

2.3  Previous  Work 

The  finite  element -alternating  method  has  its  basis  in  the  Schwartz- 
Newman  method  as  reported  by  Kantorovich  and  Krylov  (36).  Similar 
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alternating  methods  have  been  widely  used  in  the  solution  of  stress 
intensity  factor  problems  (13,  37,  38,  39). 

Prior  to  Browning  (3),  the  modeling  of  Solution  1 (Figure  1.3)  was 
formed  by  overlapping  semi -inf ini te  solids  which  greatly  complicated 
the  alternating  routine  and  restricted  the  finite  body  to  a flat  plate. 
Brownings'  introduction  of  a three-dimensional  finite  element  method 
as  Solution  1 not  only  simplified  the  algorithm  but  allowed  modeling 
of  more  complex  bodies  such  as  plates  with  fastener  holes.  Wilson's 
(40)  SSAP3  finite  element  program  with  20  node  isoparametric  elements 
was  extensively  modified  by  Browning  for  this  purpose. 

The  problem  of  an  elliptical  crack  embedded  in  an  infinite  solid 
(Figure  1.4),  which  is  used  as  Solution  2 in  the  finite  element- 
alternating  method,  was  first  considered  by  Green  and  Sneddon  (41). 

They  solved  for  stresses  when  the  crack  was  subjected  to  a uniform 
normal  pressure.  Shah  and  Kobayashi  (22)  extended  the  work  of  Green 
and  Sneddon  to  the  case  of  arbitrary  normal  pressure  on  the  crack  where 
this  pressure  Is  expressible  as  a ten  term  polynominal  in  x and  y; 

p(x,y,o)  = I I A xV  (2.3.1) 

m=o  n=o 

where  m + n < 3. 

The  ten  term  limitation  is  due  to  the  extensive  and  tedious  work 
required  in  the  calculation  of  stresses.  Smith  and  Sorenson  (39) 
programned  this  solution  for  three  stresses  and  six  terms  in  the  pres- 
sure polynominal.  Ganong  (4)  extended  the  solution  to  ten  terms  and 
all  six  stresses. 
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Ganong  also  performed  some  auxiliary  studies  to  understand  the 
character  of  the  finite  element-alternating  method.  First,  he  noted 
that  Solution  1 predicted  normal  stresses  only  on  a portion  of  the 
crack  plane  where  the  full  elliptical  crack  and  the  finite  body  over- 
lapped. Since  Solution  2 requires  definition  of  the  distribution  over 
the  entire  elliptical  crack  surface,  the  question  of  how  to  describe 
the  pressure  on  the  portion  of  the  crack  lying  "outside"  the  finite 
body  had  to  be  answered.  Ganong  chose  the  edge  crack  problem  as  a test 
case.  Five  different  valid  pressure  distributions  were  applied  to  the 
"outside"  section  of  the  crack.  Results  for  all  cases  agreed  to  within 
1%  but  convergence  was  quite  different.  To  study  the  accuracy  of  the 
polynominal  fit  to  these  crack  pressures,  he  ran  the  same  test  cases 
for  19  iterations.  Based  on  the  best  accuracy  and  most  rapid  con- 
vergence, a description  of  the  crack  pressure  "outside"  was  chosen  and 
labeled  "linear  decay."  This  produced  a pressure  distribution  "outside" 
which  linearly  decays  toward  zero  based  upon  the  pressure  along  the 
crack-finite  body  boundaries  and  the  crack  dimensions. 

Ganong  also  established  the  convergence  and  accuracy  of  the  finite 
element  portion  of  the  solution.  He  applied  uniaxial  loading  to  meshes 
of  16,  45,  96  and  112  elements  and  compared  tangential  stresses  on  the 
plane  z=0  to  those  due  to  the  analytic  solution  of  Sternberg  and 
Sadowsky.  This  study  showed  that  finer  divisions  of  the  plate  improved 
accuracy  with  the  finest  division,  112  elements,  having  a maximum 
average  deviation  of  2.5%. 

2.4  Accuracy  and  Sensitivity  Test  of  the  Method 


A three-dimensional  approximation  to  the  two-dimensional  through- 
crack  problem  (Figure  1.1)  was  made  using  two  different  finite  element 
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mesh  configurations  in  Solution  1 of  the  finite  element-alternating 
method.  The  purpose  of  this  comparison  was  to  test  both  the  accuracy 
of  the  method  and  sensitivity  to  the  number  of  finite  elements  covered 
by  the  crack. 

The  first  mesh  configuration  will  be  discussed  later  in  this 
chapter  and  is  designed  to  solve  the  problem  of  a single  quarter- 
elliptical  crack  emanating  from  a fastener  hole.  It  is  pictured  in 
Figure  2.4.  The  other  mesh  (Figure  4.1)  is  as  described  in  Chapter  4 
for  use  in  the  solution  of  the  problem  of  an  embedded  crack  adjacent  to 
a fastener  hole.  This  mesh  is  an  exact  division  of  the  single  crack 
mesh  in  the  vicinity  of  the  crack.  The  result  is  twice  the  number  of 
elements  through  the  plate  thickness  and  across  the  plate  width. 

The  through-crack  was  approximated  by  a very  slender  part- 
elliptical  crack  with  a major  to  minor  axis  ratio  of  2.5  or  greater  as 
shown  in  Figure  2.3.  The  crack  center  is  located  at  the  hole  edge  in 
the  center  of  the  plate. 

In  the  finite  element  portion  of  the  solution  method  presented  in 
this  chapter,  the  order  of  the  shape  functions  determines  the  order  of 
the  displacement  field  modeled.  The  shape  functions  for  the  20-node 
isoparametric  element  are  quadratic  in  any  one  local  coordinate  direc- 
tion. Application  of  the  strain-displacement  equation  results  in  a 
linear  description  of  strain.  A constitutive  law  links  stresses  with 
strains  and  preserves  the  linear  order.  The  concern  here  is  that  the 
finite  element  method  approximates  a non-linear  stress  field  with  a 
linear  one  over  a crack  which  is  confined  to  the  boundaries  of  a small 
number  of  finite  elements  resulting  in  excessive  stress  intensity 
factor  errors.  These  comments  apply  to  the  order  of  displacements. 
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strains  and  stresses  prior  to  mapping  into  another  coordinate  system 
which  may  or  may  not  preserve  this  order.  However,  in  this  work  all 
mid-edge  nodal  points  are  located  midway  between  the  corner  nodes  which 
maintains  linearity  in  stresses  and  strains  during  transformations. 

A major  portion  of  this  possible  difficulty  is  avoided  as  follows. 
Instead  of  initially  solving  for  stresses  at  the  crack  location  in  the 
unflawed  finite  body  using  the  finite  element  method  (Step  a.  Figure 
2.2)  an  analytic  solution  provides  this  information.  Stress  gradients 
over  the  crack  from  later  iterations  are  more  linear  and  thus  the 
finite  element  approximation  is  better.  Furthermore,  these  iterations 
contribute  only  about  25%  of  the  final  answer. 

Results  of  the  accuracy  and  sensitivity  tests  are  listed  in  Table 
2.1  and  depicted  in  Figure  2.2.  Comparisons  are  made  with  the  through 
crack  results  of  Tweed  and  Rooke  (6)  for  L/R=1.5,  1.0,  and  0.5.  In  the 
finite  element-alternating  method,  as  L/R  decreases  the  total  number  of 
elements  covered  by  the  crack  decreases.  It  is  convenient  to  view  this 
as  a reduction  in  layers  of  elements  in  the  width  direction  but  not 
through  the  thickness. 

Inspection  of  Table  2.1  shows  that  deviations  from  the  Tweed  and 
Rooke  results  do  not  exceed  3.2%.  The  maximum  difference  occurs  when 
the  crack  is  positioned  over  only  one  layer  of  elements  across  the 
width  indicating  that  one  element  coverage  in  the  width  direction  gives 
sufficient  accuracy. 

The  coverage  of  layers  of  elements  through  the  plate  thickness  is 
always  four  for  the  more  coarse  single  crack  mesh  and  eight  for  the 
embedded  crack  mesh,  regardless  of  L/R.  The  stress  gradient  through 
the  thickness  Is  more  shallow  than  across  the  plate  width.  The  finite 
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elment  solution  can  produce  a better  approximation  to  this  shallow 
gradient.  If  one  element  coverage  in  the  width  direction  is  suffi- 
cient, surely  one  element  in  the  thickness  direction  is  adequate  for 
reasonable  results.  Therefore,  very  small  cracks  bounded  by  only  one 
finite  element  can  be  used  with  the  assurance  that  stress  intensity 
factors  may  be  only  slightly  less  accurate  than  those  for  larger  cracks. 

Results  for  several  small  cracks  can  be  found  in  this  work.  The 
next  section  shows  a comparison  with  a fatigue  experiment  by  another 
author.  One  crack  considered  has  a/c=l .4  and  a/t=0.3.  Appendix  C 
presents  results  for  a "benchmark"  problem  with  a/c=2.0  and  a/t=0.2. 

Both  of  these  cracks  lie  essentially  within  the  boundaries  of  one 
finite  element  and  results  are  reasonable  with  no  signs  of  ill  behavior 
in  the  solution. 

2.5  Theoretical  Solution  Results  for  Single  Quarter-Elliptical  Cracks 

Emanat ing  from  an  Open  Fastener  Hole 

Figure  2.1a  shows  a single  quarter-elliptical  corner  crack 
emanating  from  an  open  fastener  hole  in  a finite  thickness  plate.  The 
crack  center  lies  at  the  intersection  of  the  hole  and  front  surface. 

The  plate  is  subjected  to  uniaxial  tensile  loading  in  a direction  per- 
pendicular to  the  crack  plane. 

The  finite  element  mesh  to  model  Solution  1 of  the  finite  element- 
alternating  method  for  this  crack  configuration  is  shown  in  Figure  2.4 
for  a hole  diameter  to  plate  thickness  ratio,  d/t,  of  1.0  and  in  Figure 
2.5  for  d/t=0.5.  The  mesh  is  semicircular  with  112  elements  and  z=0  is 
a plane  of  symmetry.  There  are  four  layers  of  elements  through  the 
thickness  and  seven  divisions  across  the  plate  half-width.  Planes 


21 


inclined  at  10°,  20°,  and  90°  from  the  z=0  plane  further  divide  the 
plate.  Finer  element  divisions  exist  on  the  side  of  the  hole  where  the 
crack  is  to  be  located  to  provide  improved  definition  of  stress  gradi- 
ents there. 

Ganong  (4)  used  the  finite  element-alternating  method  to  find  mode- 
one  stress  intensity  factors  for  a series  of  single  quarter-elliptical 
cracks  emanating  from  an  open  fastener  hole.  In  the  interest  of  com- 
pleteness, his  results  are  presented  in  Figures  2.6  through  2.9. 

Several  new  results  for  cases  not  considered  by  Ganong  has  been  added 
and  will  be  discussed  later  in  this  chapter. 

The  calculations  performed  for  numerous  single  quarter-elliptical 
corner  cracks  emanating  from  an  open  fastener  hole  are  listed  in 
Appendix  A.  Figures  2.$  through  2.11  are  plots  of  magnification  factor 
versus  parametric  angle  from  the  front  surface  where  the  magnification 
factor  is  defined  as 
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(2.5.1) 


Kj  is  the  mode-one  stress  intensity  factor,  a is  the  crack  dimension 
along  the  hole  edge,  and  a is  the  remote  uniaxial  tensile  stress.  The 
parametric  angle,  9,  is  defined  by  the  parametric  equations  for  an 
ellipse 

X = A sin  e x = A cos  o 

or  (2.5.2) 

y = B cos  0 y = B sin  0 

where  A and  B are  the  semi -major  and  semi -minor  axes  of  the  ellipse. 
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and  X and  y are  coordinates  of  a point  on  the  periphery  of  the  ellipse. 
All  finite  element-alternating  method  calculations  in  this  work  are  for 
a plate  having  a modulus  of  elasticity  of  30x10^  psi  and  a Poisson's 
ratio  of  0.25. 

Ganong  showed  results  for  crack  shapes  of  a/c=0.99  and  1.01  which 
are  omitted  here  for  reasons  to  follow.  These  crack  shapes  were  used 
to  approximate  a circular  crack  with  a near-circular  elliptical  crack 
for  comparison  with  results  of  other  authors.  Due  to  the  limiting 
forms  of  stresses  in  the  solution  of  the  elliptical  crack  problem  (Solu- 
tion 2),  exact  circles  cannot  be  duplicated.  Ganong  reasoned  that  since 
results  for  a/c=0.99  and  a/c=1.01  were  very  close,  both  were  accurate. 
Unfortunately,  further  study  shows  that  both  cases  are  subject  to 
excessive  errors.  The  difficulties  lie  in  the  formation  of  the  matrix 
which  yields  stress  function  coefficients  in  Solution  2 of  the  finite 
element-alternating  method.  Inspection  of  the  equations  describing 
terms  in  this  matrix  shows  definite  tendencies  for  large  round-off 
errors  in  the  numerical  calculations  when  the  major  and  minor  axes  of 
the  elliptical  crack  are  nearly  equal.  Shah  and  Kobayashi  (22)  alluded 
to  this  round-off  error  problem  when  the  circular  crack  shape  was 
approximated  too  closely  by  a near-circular  elliptical  crack.  As  a 
result,  they  simulated  a penny-shaped  crack  with  an  elliptical  crack 
of  shape  a/c=0.982.  In  the  finite  element-alternating  method,  the 
improper  formation  of  the  matrix  results  in  incorrect  values  of  the 
stress  function  coefficients.  This,  in  turn,  causes  the  calculation  of 
Incorrect  stresses  on  the  surface  of  the  finite  body  by  Solution  2. 

When  these  stresses  are  reversed  to  free  the  surfaces  and  applied  to 
the  finite  body  to  find  the  new  crack  plane  stress  distribution,  the 
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body  has  a tendency  to  attempt  a large  rigid  body  translation.  Because 
one  nodal  point  is  always  fixed  in  x and  y,  (0,  0,  6)  in  this  case,  the 
body  cannot  translate  but  bends  instead.  The  result  is  large  errors  in 
the  description  of  the  crack  plane  stresses  and  incorrect  stress  in- 
tensity factors.  As  an  example  of  the  magnitude  of  these  errors,  if 
Ganong's  results  for  a single  quarter-elliptical  crack  of  a/c=l .01  and 
a/t=0.5  are  added  to  the  parametric  plot  in  Figure  6.24,  a 30%  error  in 
the  value  of  the  stress  intensity  factor  at  the  front  surface  is  seen. 
Calculations  have  been  successfully  carried  out  for  crack  shapes  of 
a/c=0.90  and  1.10.  No  calculations  have  been  completed  in  the  range 
0.90<a/c<0.99,  yet  it  is  thought  that  due  to  Shah  and  Kobayashi's 
observation,  the  shapes  to  avoid  are  probably  0.98<a/c<l .02. 

Ganong  (4)  identified  six  effects  which  control  the  variation  of 
magnification  factor  around  the  crack  periphery.  He  interpreted  the 
results  in  Figures  2.6  to  2.9  in  terms  of  these  six  effects.  Except 
for  some  general  observations,  his  detailed  comments  will  not  be  re- 
peated. However,  the  six  effects  will  be  described  here  because  they 
will  be  useful  later. 

1.  Hole  Stress  Concentration  Effect.  Stress  intensity  factors 
are  sensitive  to  the  local  stress  field  in  the  vicinity  of 
the  crack.  Cracks  having  large  c/r  (Figure  2.1a),  where  r 
is  the  hole  radius,  extend  into  regions  of  lower  stresses 
away  from  the  hole  and  would  be  expected  to  have  lower  stress 
intensity  factors.  The  reverse  is  true  when  low  c/r  places 
the  entire  crack  In  an  area  of  elevated  stress  and  stress 
Intensity  factors  rise. 
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2.  Front  Surface  Effect.  Cracks  intersecting  a free  surface  can 
open  more  than  cracks  embedded  in  a solid  under  identical 
loading.  Increased  crack  opening  is  reflected  in  higher 
stress  intensity  factors. 

3.  Back  Surface  Effect.  As  the  ratio  a/t  increases  and  the 
crack  approaches  the  back  surface,  stress  intensity  factors 
increase.  This  is  due  to  the  decreased  rigidity  of  the 
ligament  between  the  crack  and  back  surface. 

4.  Hole  Surface  Effect.  For  small  cracks,  the  hole  surface  acts 
like  a free  surface.  This  allows  the  crack  to  open  and  conse- 
quently increases  stress  intensity  factors.  For  larger  cracks, 
the  hole  acts  to  restrict  crack  opening  and  stress  intensity 

is  lowered. 

5.  Crack  Shape  Effect.  The  largest  stress  intensity  factor 
for  an  elliptical  crack  under  uniform  loading  occurs  at  the 
ends  of  the  minor  axis.  The  crack  shape  effect  decreases 
as  an  elliptical  crack  approaches  a circular  shape. 

6.  Through-Thickness  Stress  Variation  Effect.  For  a hole  diameter 
to  plate  thickness  ratio  of  1.0,  the  stress  concentration  at 
the  hole  edge  varies  from  2.7  at  the  plate  surfaces  to  3.1 

at  the  center  of  the  plate.  A similar  variation  occurs 
through  the  thickness  at  locations  removed  from  the  hole  and 
results  in  like  trends  in  stress  intensity  factors. 

In  the  present  study.  It  is  not  possible  to  separate  the  six 
effects  and  show  their  Importance  Independent  of  one  another.  To  do 
so  would  require  corresponding  results  for  a problem  In  which  one  of 
the  effects  was  not  Important.  Attempts  to  accomplish  this  usually 
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result  in  two  or  more  effects  changing  simultaneously  rather  than  just 
one.  Only  the  shape  effect  can  be  so  isolated  at  present  as  suggested 
by  Kobayashi  (42).  He  normalizes  the  magnification  factor  (Equation 
2.5.1)  by  multiplying  by  N,  where 

N = aV(A^cos20  + B^sin^e)^  , (2.5.3) 

A and  B are  the  semi-major  and  semi-minor  axes  of  the  elliptical  crack 
and  0 is  the  parametric  angle  from  the  minor  axis. 

Figures  2.6  through  2,9  are  plots  at  constant  a/t  for  varying 
crack  shapes,  a/c.  At  the  front  surface  intersection,  magnification 
factors  change  no  more  than  about  8%  over  a wide  range  of  crack  shapes 
due  to  a balance  maintained  between  a number  of  effects.  Increasing 
di/c  means  decreased  c/r  which  raises  stress  intensity  factors  due  to 
the  hole  surface  effect.  Decreasing  c/r  also  places  the  crack  border 
near  the  front  surface  in  a region  of  higher  stresses  which  leads  to 
higher  stress  intensities.  For  cracks  with  a/c<3.0,  the  increase  in 
magnification  factor  as  the  front  surface  is  approached  is  due  to  a 
combination  of  front  surface  and  crack  shape  effects. 

Ganong  (4)  remarked  that  magnification  factors  were  always  highest 
at  the  hole  intersection  due  to  the  dominant  character  of  the  hole 
stress  concentration  effect.  While  this  may  be  true  for  single  quarter- 
elliptical  cracks  of  a/c<2.0  as  considered  by  Ganong,  it  is  not  the 
case  for  a/c>2.0.  Figures  z.6  and  2.8  show  magnification  factors  for 
a/c=3.0  and  a/t  of  0.9  and  0.5.  In  both  cases,  the  location  of  maximum 
stress  intensity  is  shifted  to  an  intermediate  point  along  the  crack 
periphery  and  nearer  the  front  surface.  This  is  again  due  to  a combina- 
tion of  several  effects.  As  a/c  increases  to  3.0,  c/r  decreases  and 
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locations  on  the  crack  border  near  the  front  surface  He  in  an  area  of 
higher  stresses.  Also,  the  crack  shape  effect  causes  stress  intensity 
factors  to  decrease  at  the  hole  intersection. 

The  variation  in  magnification  factor  along  the  crack  periphery 
for  cracks  with  a/c<1.0  is  wider  than  for  a/c>l .0  as  can  be  seen  in 
Figures  2.8,  2.9  and  2.10.  Since  the  crack  minor  axis  is  along  the 
front  surface  for  a/c<l ,0,  the  crack  shape  effect  does  not  reduce  the 
elevated  stress  intensity  at  the  hole  intersection  as  happens  with 
higher  a/c  values.  The  shape  effect,  hole  surface  effect  and  hole 
stress  concentration  effect  combine  to  lower  stress  intensity  factors 
at  the  front  surface  for  cracks  which  are  long  in  that  direction. 

Figures  2.9  and  2.10  show  this  to  be  strong  enough  to  overcome  the 
effect  of  the  front  surface. 

Figure  2.11  shows  magnification  factor  versus  parametric  angle  for 
the  same  crack  with  hole  diameter  to  plate  thickness  ratios  of  1.0  and 
0.5.  This  reduction  was  accomplished  by  halving  the  hole  diameter 
while  leaving  all  other  plate  and  crack  dimensions  unchanged.  Thus, 
the  only  effects  on  the  magnification  factor  should  be  those  of  the 
hole  surface  and  the  hole  stress  concentration.  As  expected,  the  hole 
stress  concentration  effect  dominates,  resulting  in  a lowering  of  mag- 
nification factor  everywhere  as  the  increased  c/r,  due  solely  to  reduced 
hole  radius  , places  the  crack  in  a region  of  lower  stresses.  The 
effect  of  the  hole  surface  can  be  seen  in  the  wider  separation  of  the 
two  curves  at  the  hole  intersection.  The  reduction  in  c/r  increases 
the  rigidity  of  the  hole  with  respect  to  the  unchanged  crack  size. 

This  increased  rigidity  is  reflected  in  reduced  crack  opening  and 
lowered  stress  Intensity  factors  for  d/t=0.5. 
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I * It  is  perhaps  important  to  recognize  that  stress  intensity  factors 

j presented  in  this  work  at  the  hole  and  surface  intersections  may  be 

somewhat  in  error.  Hartranft  and  Sih  (43)  studied  the  semi-circular 
I surface  flaw  in  a semi-infinite  solid  in  great  detail  and  discovered 

that  stress  intensity  factors  on  the  crack  periphery  near  the  surface 
^ dropped  sharply  toward  zero.  Sih  (44)  explains  this  result  with  the 

j I observation  that  the  state  of  stress  is  very  complex  near  a surface. 

I At  all  other  locations,  Hartranft  and  Sih's  results  compare  closely 

I with  those  of  Smith  (45),  This  decrease  in  stress  intensity  factor 

near  the  surface  intersections  was  not  seen  in  the  present  work,  due 
I probably  to  the  nature  of  the  approximations  used.  However,  experiments 

I conducted  by  the  author  show  cracks  to  propagate  at  points  of  highest 

^ stress  intensity  factor,  not  lowest.  A small  shift  in  this  location, 

while  theoretically  important,  has  little  practical  significance. 

2.6  Comparison  of  Theoretical  Results  with  Two  Independent  Experiments 
Snow  (20)  reports  a series  of  seven  tests  resulting  in  predicted 
stress  intensity  factors  for  some  single  quarter-elliptical  cracks  at 
open  fastener  holes  (Figure  2.1a)„  In  his  work,  a starter  crack  was 

f 

introduced  at  the  intersection  of  a hole  and  a surface  in  a plate  of 
~ polymethylmethacrylate  (PMMA).  The  plate  was  cycled  under  low  load 

levels  and  the  crack  grew  by  fatigue.  At  periodic  intervals,  35  mm 
> photographs  were  taken  perpendicular  to  the  crack  plane.  Measurements 

r,  from  the  photographs  gave  crack  length  versus  number  of  cycles  data. 

The  Paris  fatigue  relationship, 

^*C(AK)'"  (2.6.1) 

1 was  assumed  to  apply  where  C and  m are  material  dependent  constants. 


11 


28 


I 


I 

( 

( 

I 


' r 

I 

i 

\ 


aK  is  the  change  in  stress  intensity  factor  during  cycling,  and  da/dn 
is  the  crack  growth  rate.  C and  m were  determined  from  a single 
through-crack  baseline  test  (Figure  1.1)  where  da/dn  was  measured  and 
aK  was  known.  Crack  length  versus  cycles  data  from  the  fastener  hole 
tests  required  careful  differentiation  to  yield  crack  growth  rates. 
Knowledge  of  da/dn,  C,  and  m permits  solution  of  the  Paris  relationship 
for  stress  intensities. 

The  most  complete  results  shown  by  Snow  are  for  an  average  a/c  of 
1.4  and  a range  of  a/t  from  0.23  to  0.84.  Magnification  factors  are 
presented  at  the  hole  and  surface  intersections  only.  The  average  hole 
diameter  to  plate  thickness  ratio,  d/t,  for  the  seven  tests  was  1.05. 
This  close  agreement  in  plate  geometries  between  Snow's  experiment  and 
the  finite  element-alternating  method  model  permits  a valid  comparison 
of  results  from  the  two  studies. 

To  form  this  comparison,  stress  intensity  factors  were  calculated 
by  the  finite  element-alternating  method  for  single  quarter-elliptical 
cracks  emanating  from  an  open  fastener  hole  where  a/c=1.4,  d/t=1.0,  and 
a/t  ranged  from  0.3  to  0.8.  Figures  2.12  and  2.13  are  plots  of  magnifi- 
cation factor  versus  a/t  for  the  experiments  of  Snow  and  the  theoretical 
results  of  the  finite  element-alternating  method.  Agreement  at  the  hole 
intersection  is  within  2%.  Agreement  at  the  front  surface  intersection 
is  not  as  close  but  trends  with  increasing  a/t  are  identical.  Here, 
the  finite  element-alternating  method  predicts  magnification  factors 
which  average  13%  below  Snow's  results. 

Figures  2.14  and  2.15  show  a similar  comparison  for  an  average  a/c 
of  1.5.  Agreement  between  Snow's  results  and  those  of  the  finite 
element-alternating  method  is  with  3%  at  both  the  hole  and  surface 
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intersections  for  a/t=0.5  and  trends  with  increasing  a/t  are  in  close 
accord. 

The  Poisson's  ratio  of  PMMA  is  0.35  while  all  finite  element- 
alternating  method  calculations  are  for  a value  of  0.25.  Smith  (46) 
noted  a 3%  increase  in  back  surface  magnification  factors  when  Pois- 
son's ratio  was  increased  from  0.25  to  0.39  in  his  theoretical  study 
of  semi-elliptical  surface  cracks.  If  all  theoretical  results  shown 
in  Figures  2.12  through  2.15  are  increased  by  this  amount,  the  max- 
imum disagreement  is  reduced  to  lOili. 

McGowan  and  Smith  (19)  used  stress-freezing  photoelasticity  to 
experimentally  determine  stress  intensity  factors.  Their  results  were 
limited  to  the  points  of  intersection  of  the  crack  with  the  hole  and 
front  surface.  Table  2.2  shows  a comparison  of  magnification  factors 
from  McGowan  and  Smith's  experiments  with  results  from  the  finite 
element-alternating  method.  Only  three  of  the  eight  tests  have  plate 
and  crack  geometries  close  to  those  of  existing  finite  element- 
alternating  method  results.  The  striking  feature  of  most  of  the 
experimental  tests  is  a higher  magnification  factor  at  the  front 
surface  than  at  the  hole.  McGowan  and  Smith  claim  this  is  due  to  load 
transfer  to  the  side  of  the  hole  opposite  the  crack  which  tends  to 
lower  the  magnification  factor  at  the  hole  intersection.  Test  4 
shows  this  lowered  value  for  a crack  with  a/c=1.10  and  a/t=0.48. 

Figure  6.25  is  a parametric  plot  of  stress  intensity  factor  versus 
various  values  of  a/c  and  a/t  from  finite  element-alternating  method 
calculations  to  explain  slow  growth  in  the  experiments  discussed  in 
Chapter  6.  This  plot  shows  that  the  change  in  location  of  maximum 
magnification  factor  for  cracks  with  a/t=0.5  should  not  occur  until 
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Surface 

Hole 

Test  4 

0.48 

1.10 

0.88 

0.90 

0.89 

Test  6 

0.46 

1 .55 

0.94 

0.98 

0.94 

Test  7 

0.74 

1.98 

0.95 

1 .08 

0.66 

Finite  Element- 
Alternating  Method 
(equivalent  to) 

Test  4 

0.50 

1.10 

1.00 

1 .03 

1 .36 

Test  6 

0.50 

1 .50 

1 .00 

1 .04 

1 .24 

Test  7 

0.75 

2.00 

1.00 

0.94 

1.06 
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® a/c  is  larger  than  2.0  This  suggests  that  McGowan  and  Smith  may  be 

I observing  a valid  trend  but  that  onset  of  the  shift  in  location  of 

maximum  magnification  factor  is  being  seen  at  crack  shapes  that  are 
I much  too  small. 

, Static  fracture  tests  of  plates  having  cracks  of  various  shapes 

^ and  locations  relative  to  fastener  holes  were  performed  by  the  author, 

y Comparison  of  the  present  theory  with  a large  number  of  these  exper- 

iments is  presented  in  Chapter  6. 
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CHAPTER  3 


DOUBLE  SURFACE  CRACKS  EMANATING  FROM  AN  OPEN  FASTENER  HOLE 

3. 1 Problem  Statement 

Figure  2.1b  shows  double  quarter-elliptical  surface  cracks 
emanating  from  an  open  fastener  hole  in  a finite-thickness  plate.  The 
crack  centers  lie  at  the  intersection  of  the  hole  and  the  front  surface 
The  plate  is  subject  to  uniaxial  tensile  loading  in  a direction  per- 
pendicular to  the  plane  containing  both  cracks.  Regardless  of  the 
orientation  of  the  crack  major  and  minor  axes,  "a"  refers  to  that  crack 
dimension  along  the  edge  of  the  hole  while  "c"  is  along  the  front  sur- 
face of  the  plate.  Both  cracks  are  of  the  same  shape  and  depth.  This 
double  crack  configuration  is  seen  in  fatigue  failures  when,  following 
initiation  of  a single  crack,  a second  crack  propagates  on  the  side  of 
the  hole  opposite  the  first.  Mode  one  stress  intensity  factors  are  to 
be  calculated  around  the  peripheries  of  both  cracks. 

3.2  Finite  Element  Mesh  Configuration 

The  finite  element  mesh  used  in  the  solution  of  the  double  crack 
problem  is  quarter-circular  with  108  elements  (Figure  3.1).  This  mesh 
is  arranged  similar  to  those  for  the  single  and  embedded  crack  config- 
urations. There  are  four  layers  of  elements  through  the  thickness  and 
nine  divisions  across  the  plate  half-width.  Planes  inclined  at  10°  and 
20°  from  the  plane  z-O  further  divide  the  plate.  Finer  element  divi- 
sions exist  near  the  hole  to  provide  improved  definition  of  the  stress 
gradients  there.  The  planes  z=0  and  y=0  are  assumed  planes  of  symmetry 
and  are  rollered  to  prevent  displacements  in  the  z and  y directions. 
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respectively.  The  mesh  is  pictured  in  Figure  3.1  for  a hole  diameter 
to  plate  thickness  ratio,  d/t,  of  1.0  and  in  Figure  3.2  for  d/t=0.5. 

3.3  Changes  to  the  Finite  Element-Alternating  Method  for  the  Double 

Crack  Problem 

The  finite  element-alternating  method  as  applied  to  the  problem  of 
single  quarter-ell iptical  cracks  emanating  from  a fastener  hole  uses  as 
one  solution  that  of  a single  elliptical  crack  in  an  infinite  medium 
subjected  to  an  arbitrary  pressure  distribution.  No  solution  is  cur- 
rently available  for  the  companion  problem  of  two  cracks.  The  diffi- 
culties caused  by  use  of  the  single  crack  solution  coupled  with  the 
finite  element  mesh  of  Figures  3.1  or  3.2  for  the  double  crack  problem 
stem  from  a mismatch  in  planes  of  symmetry.  While  symmetry  about  the 
y=0  plane  is  forced  upon  the  finite  element  solution,  the  crack  solu- 
tion has  no  such  symmetry  about  this  plane.  The  only  symmetries 
recognized  by  the  crack  solution  are  those  about  planes  passing  through 
the  crack  center,  none  of  which  coincide  with  the  plane  y=0  because  the 
crack  center  is  offset  from  the  hole  center.  Therefore,  changes  to  the 
finite  element-alternating  method  were  necessary  to  allow  for  the 
presence  of  the  second  crack  and  to  produce  a match  in  the  symmetries. 

The  following  changes  were  made  to  the  solution  method  for  the 
double  crack  problem  and  represent  one  iteration.  The  crack  referred 
to  as  the  "right"  crack  has  its  center  at  the  intersection  of  the  fas- 
tener hole  and  front  surface  on  the  side  of  the  plate  actually  modeled 
by  the  finite  element  mesh.  The  "left"  crack  is  a mirror  image  of  the 
"right"  crack  about  the  y=0  plane.  Between  Steps  a and  b.  Figure  2.2, 
Solution  2 of  the  finite  element-alternating  method  calculates  residual 
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stresses  at  the  location  of  the  left  crack  due  to  the  freeing  of  resid- 
ual stress  on  the  right  crack.  These  residuals  are  combined  and  applied 
to  the  right  crack  (Step  b.  Figure  2.2).  Solution  2 then  solves  for 
stresses  at  locations  on  the  finite  body  and  mirror  image  locations  on 
a fictitious  quarter  plate  mirrored  about  the  y=0  plane  (Figure  3.1). 
Combined  stresses  from  actual  and  mirror  locations  are  residuals  to  be 
removed  from  the  finite  body  surfaces  (Step  c.  Figure  2.2). 

Some  care  must  be  exercised  in  the  use  of  this  double  crack  al- 
gorithm. If  the  major  or  minor  axis  of  one  crack  touches  the  center  of 
the  second  crack,  the  elliptical  crack  solution  is  asked  to  calculate 
singular  stresses  on  the  crack  boundary.  The  obvious  correction  is  to 
avoid  crack  dimensions  which  create  this  problem.  At  first  glance,  the 
overlapping  of  the  two  cracks  also  raises  questions  of  solution  valid- 
ity. Crack  overlap  causes  the  elliptical  crack  solution  to  calculate 
stresses  on  the  plane  z=0  inside  the  crack  which  in  itself  causes  no 
problems.  Results  presented  in  this  thesis  are  for  cracks  which  over- 
lap only  in  regions  outside  the  boundaries  of  the  finite  body.  Since 
the  stresses  in  these  regions  are  linearly  decayed  prior  to  the  poly- 
nomial fitting  procedure  which  describes  the  crack  plane  stress 
distribution,  any  possible  problems  due  to  crack  overlap  are  eliminated. 

3.4  Three  Dimensional  Analysis  of  a Two-Dimensional  Test  Problem 

Bowie  (5)  solved  the  two-dimensional  problem  of  double  through- 
cracks  emanating  from  a hole  in  a plate  subjected  to  uniaxial  loading 
(Figure  1.2).  His  well  respected  solution  employed  complex  variable 
techniques  to  calculate  stress  intensity  factors.  To  form  a comparison 
using  the  present  method,  double  through-cracks  were  approximated  by 
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two  very  slender  part-el  1 iptical  cracks  with  major  to  minor  axis  ratios 
of  2.5  or  greater.  This  configuration  is  pictured  in  Figure  3.3  along 
with  results  from  the  finite  element-alternating  method  and  Bowie's 
solution.  Calculations  were  performed  for  L/R  of  1.5,  1.0,  and  0.5. 
Results  at  the  center  of  the  plate  are  presented  to  minimize  the  effects 
of  the  surfaces  and  thus  more  closely  model  a plane  strain  situation. 

The  three-dimensional  finite  element-alternating  method  analysis  of  this 
two-dimensional  problem  shows  agreement  to  within  5%  of  the  Bowie 
results. 

3.5  Solution  Results 

A list  of  calculations  performed  for  various  double  cracks 
emanating  from  open  holes  is  shown  in  Appendix  A.  Comparitive  plots  of 
magnification  factor  versus  parametric  angle  from  the  front  surface  are 
presented  in  Figures  3.4  through  3.9.  The  magnification  factor  is 
defined  as  before. 


Crack  shapes  and  depths  used  in  the  double  crack  calculations  were 
identical  to  those  used  by  Ganong  (4)  for  single  cracks.  This  was 
done  with  forethought  so  that  effects  of  the  addition  of  a second  crack 
could  be  analyzed.  The  results  for  the  corresponding  single  cracks  are 
shown  and  discussed  in  Chapter  2 and  ^orm  a basis  for  comparison  of 
double  crack  results.  Both  this  author  anc<  Ganong  analyzed  the  single 
crack  results  with  respect  to  the  six  effects  mentioned  in  Chapter  2. 
These  comments  apply  equally  well  to  the  double  crack  results  with  the 
possible  exception  of  the  hole  surface  effect. 


36 


I 

I 

I 

t 

r 

I 

*- 

( 

I 

I 

I 

f 

\ 

i 

f 

I 

I 

f 

r» 


For  small  single  cracks,  the  hole  surface  acts  like  a free  surface 
which  tends  to  elevate  stress  intensity  factors.  Bowie  (5)  noted  for 
through-cracks  that  when  the  crack  length  to  hole  radius  ratio  is  small, 
stress  intensity  factors  for  one  and  two  cracks  are  approximately  the 
same.  As  cracks  get  larger,  however,  the  effect  of  the  hole  is  quite 
different  for  single  and  double  cracks.  For  larger  single  cracks,  hole 
rigidity  restricts  crack  opening  which  lowers  stress  intensity  factors. 
The  addition  of  a second  large  crack  has  the  opposite  effect  since  it 
reduces  hole  rigidity  and  allows  both  cracks  to  open  more.  This  pro- 
duces elevated  stress  intensity  factors  with  respect  to  the  single  crack 
case.  All  double  crack  calculations  presented  here  exhibit  this  trend. 

Casual  inspection  of  plots  of  magnification  factor  versus  para- 
metric angle  (Figures  3.4  through  3.9)  shows  little  difference  from  the 
single  crack  results  (Figures  2.6  through  2.10).  While  differences  do 
appear  subtle,  they  can  be  as  large  as  8%.  In  all  cases,  the  maximum 
elevation  in  stress  intensity  factors  with  respect  to  the  single  crack 
occurred  not  at  crack  intersections  with  surfaces  but  rather  at  para- 
metric angles  of  from  30  to  50  degrees.  At  these  locations,  stress 
intensity  factors  for  double  cracks  were  4-8%  higher  than  for  the  cor- 
responding single  cracks.  The  amount  of  increase  appears  to  be  a 
function  of  the  crack  length  along  the  front  surface  versus  hole  radius 
ratio,  c/r,  regardless  of  crack  shape  or  depth.  As  c/r  increases  from 
0.5  to  1.0,  the  difference  in  stress  intensity  factors  for  double 
versus  single  cracks  increases  uniformly  from  4 to  8%  at  the  30  to  50 
degree  location.  From  c/r  of  1.0  to  1.34  the  difference  remained  at 
the  high  of  7 to  8%. 


I 


37 


I 

I 

I 

( 

V 

» 

I 

I 

I 

I 

I 

I 

I 

\ 

I 

I 

r 


Increases  in  stress  intensity  factors  at  the  hole  and  front  surface 
intersections  for  double  versus  single  cracks  were  lower  than  at  the 
30-50  degree  location.  The  amount  of  change  was  a definite  function  of 
a/t  as  might  be  expected,  since  as  a/t  increases  hole  rigidity  de- 
creases. For  a/t=0. 5, double  crack  results  show  no  increase  at  the  hole 
intersection  but  a 3-5%  increase  at  the  front  surface.  As  a/t  becomes 
larger,  stress  intensity  factors  at  the  hole  intersection  begin  to  in- 
crease until  they  are  2-5%  above  single  crack  results  while  values  at 
the  front  surface  remain  about  3-5%  higher. 

Figure  3.9  shows  magnification  factors  for  the  same  double  crack 
shape  and  hole  diameter  to  plate  thickness  ratios,  d/t,  of  1.0  and 
0.5.  As  with  the  corresponding  single  crack  problem  discussed  in 
Chapter  2,  the  only  dimensional  change  in  the  reduction  of  d/t  from 
1.0  to  0.5  was  the  hole  diameter.  Both  cracks  extend  like  distances 
along  the  hole  and  front  surfaces  yet  the  ratio  c/r  is  doubled  for  the 
d/t=0.5  case.  The  depth  ratio,  a/t,  remains  the  same.  Earlier  obser- 
vations on  double  versus  single  crack  results  linked  magnification 
factor  differences  at  the  hole  and  surface  intersections  to  a/t  and 
at  intermediate  points  to  c/r.  The  same  trend  in  reverse  seem  to  be 
occurring  here  for  changing  c/r.  As  this  ratio  increases,  magnification 
factors  decrease  at  intermediate  points,  but  do  not  change  appreciably 
near  the  hole  and  front  surfaces.  Again,  this  lends  support  to  the 
dominant  effect  of  the  hole  stress  concentration  which  overshadows 
effects  of  the  hole  surface.  Increasing  c/r  for  double  cracks  would 
be  expected  to  elevate  stress  intensity  factors  due  to  decreased  hole 
rigidity  yet  the  fact  that  the  crack  lies  in  a region  of  lowered  stress 
is  still  more  important. 
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3.6  Comparison  with  an  Estimate  of  Another  Author 

Shah  (12)  used  a Green's  function  technique  to  find  stress  in- 
tensity factors  for  two  semi-el  1 iptical  embedded  cracks  emanating  from 
a fastener  hole.  He  applied  correction  factors  for  the  surface  effects 
and  reasoned  that  for  through-cracks  (Figure  1.2) 
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(3.6.1) 


The  corresponding  elliptical  crack  relationship  was  assumed  to  be 
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where  is  the  crack  length  of  a through-crack  having  the  same  area 
as  an  elliptical  crack  of  depth  a and  length  c.  Shah's  estimates  could 
be  discussed  at  length  in  a comparison  with  solutions  by  the  finite 
element-alternating  method.  This  author  chooses  to  deal  solely  with 
the  one-to-two  crack  corrections  predicted  by  both  methods.  Shah's 
factor  actually  corrects  double  crack  results  to  find  stress  intensity 
factors  for  a single  crack  but  comments  below  will  treat  this  in  reverse 
for  comparison  with  the  present  study. 

While  Shah's  one-to-two  crack  correction  factor  tends  to  over- 
estimate the  average  differences  found  with  the  finite  element -alter- 
nating method  by  only  about  1 to  2%,  his  factor  predicts  a constant 
correction  for  all  locations  on  the  crack  front.  The  present  study 
indicates  there  can  be  as  much  as  an  8%  variation  along  the  crack 
border.  His  factor  is  also  necessarily  limited  to  a/c  < 1.0,  not  only 
because  these  are  the  shapes  he  considered  but  because  crack  orientation 
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I ^ is  important.  Consider  two  cracks  with  identical  areas  but  one  having 

^ I the  major  axis  along  the  front  surface  (a/c<1.0)  and  the  other  with  the 

! major  axis  along  the  hole  (a/c>1.0).  The  finite  element-alternating 

1 

! I method  predicts  correction  factors  differing  by  as  much  as  5%  for  these 

two  cases  while  Shah's  correction  factors  would  be  identical.  Figure 
^ 3.7  shows  a comparison  of  Shah's  estimate  with  present  results  for 

j-  double  cracks  having  a/c=0.5  and  a/t=0.25. 

Static  fracture  tests  of  plates  having  cracks  of  various  shapes 

r 

I and  locations  relative  to  fastener  holes  were  performed  by  the  author. 

Comparison  of  the  theory  with  a large  number  of  these  experiments  is 
f presented  in  Chapter  6. 
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CHAPTER  4 


EMBEDDED  CRACKS  EMANATING  FROM  AN  OPEN  FASTENER  HOLE 

4.1  Problem  Statement 

Figure  2.1c  shows  a single  half-elliptical  crack  emanating  from  an 
open  fastener  hole  in  a finite-thickness  plate.  The  crack  center  is 
located  at  the  edge  of  the  hole,  midway  through  the  plate  thickness. 

The  crack  major  axis  is  always  parallel  to  the  hole  centerline.  The 
crack  has  no  front  or  back  surface  intersections  but  some  results  are 
shown  for  cases  where  2c=t.  This  configuration  will  be  referred  to  as 
an  embedded  crack.  The  finite  width  plate  is  subject  to  uniaxial 
tensile  loading  in  a direction  perpendicular  to  the  crack  plane.  Stress 
Intensity  factors  are  to  be  calculated  along  the  crack  periphery.  While 
many  in-service  failures  have  resulted  from  corner  cracks  at  fastener 
holes,  the  embedded  crack  has  recently  gained  attention  due  to  failures 
in  aircraft  jet  engine  turbine  rotors  due  to  cracks  of  this  type. 

4.2  Fi nUe  Element  Mesh  Configuration 

The  finite  element  mesh  used  in  the  solution  of  the  embedded  crack 
problem  (Figure  4.1)  is  similar  to  that  used  by  Ganong  (4)  for  the 
single  corner  crack  emanating  from  a fastener  hole.  Both  meshes  have 
112  elements  with  the  difference  that  the  embedded  crack  mesh  has 
synmetry  about  the  x=0  plane  in  addition  to  the  z=0  plane.  Therefore, 
all  112  elements  lie  in  the  direction  of  the  positive  x-axis  and  are 
distributed  over  the  plate  half-thickness.  This  mesh  arrangement  1s 
pictured  In  Figure  4.1  for  a hole  diameter  to  plate  thickness  ratio, 
d/t,  of  1.0  and  In  Figure  4.2  for  d/t=0.5. 
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4 . 3 Solution  Results 

A list  of  calculations  performed  for  numerous  embedded  cracks 
emanating  from  an  open  fastener  hole  is  shown  in  Appendix  A.  Compara- 
tive plots  of  magnification  factor  versus  parametric  angle  from  the 
crack  minor  axis  are  presented  in  Figures  4.3  through  4.7.  The 
magnification  factor  is  defined  as 

MF  = — ^ (4.3.1) 

a/iia" 

or  for  ease  of  comparison  in  some  of  the  plots. 


o 


where  Kj,  o and  a are  as  defined  earlier. 

The  six  effects  assumed  to  alter  the  stress  intensity  factor  around 
the  crack  periphery  for  single  and  double  corner  cracks  emanating  from 
a fastener  hole,  in  general,  also  apply  to  the  embedded  crack.  These 
effects  are  listed  in  Chapter  2.  Two  important  differences  must  be 
noted,  however.  The  ends  of  the  crack  major  axis  approach  the  front  and 
back  surfaces  of  the  plate  equally.  The  effect  here  is  the  same  as  the 
back  surface  effect  for  corner  cracks  and  will  be  referred  to  simply  as 
the  surface  effect.  Also,  since  the  crack  never  intersects  the  plate 
surfaces  in  the  results  presented  here,  the  front  surface  effect  does 
not  apply  to  the  embedded  crack  configuration. 

In  all  results  shown,  the  effect  of  the  hole  stress  concentration 
is  of  great  importance  and  accounts  for  a large  share  of  the  variation 
In  magnification  factor  along  the  crack  front.  Thus,  portions  of  the 
crack  nearest  the  hole  would  be  expected  to  have  elevated  stress  Inten- 
sities while  those  removed  from  the  hole  would  have  smaller  values. 
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The  following  discussion  deals  more  with  deviations  from  this  behavior 
caused  by  other  effects  than  with  this  effect  directly. 

Figures  4.3  and  4.4  show  comparative  plots  for  five  crack  shape 
and  two  crack  depth  ratios.  Since  the  crack  depth  to  hole  radius  ratio, 
a/r,  is  constant  in  each  figure,  the  hole  surface  effect  is  deemed  not 
a factor  in  comparisons  between  results  on  the  same  figure.  Figure  4.3 
shows  the  dominant  character  of  the  crack  shape  effect.  The  surface 
and  hole  stress  concentration  effects  expected  to  occur  at  the  hole  in- 
tersections are  overshadowed  by  the  strong  shape  effect  which  tends  to 
reduce  stress  intensity  at  the  ends  of  the  crack  major  axis.  Addition- 
ally, the  magnification  factor  changes  only  slightly  at  this  location 
with  rather  large  changes  in  crack  shape.  This  is  the  result  of  a 
balance  as  a/c  increases  between: 

1.  Decreasing  surface  effect  (Kj  decreases) 

2.  Increasing  shape  effect  (Kj  increases) 

3.  Increasing  through-thickness  stress  concentration 
(Kj  increases) 

Figure  4.4  shows  trends  similar  to  those  noted  in  Figure  4.3 
except  that  shape  effects  are  now  smaller  due  to  higher  a/c  values. 

Here  the  surface  effect  is  more  important  as  seen  by  the  decrease  in 
magnification  factor  at  the  hole  intersection  as  the  crack  major  axis 
moves  away  from  the  surfaces  (increasing  a/c).  Again,  a balance  between 
the  same  effects  noted  in  the  previous  paragraph  maintains  identical 
magnification  factor  values  at  the  hole  Intersection  as  a/c  increases 
from  0.67  to  0.75. 

Figures  4.5  and  4.6  show  comparative  plots  for  two  a/c  and  four 
a/t  ratios.  Since  the  crack  shape  is  constant  in  each  figure,  the  crack 
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shape  effect  is  not  a factor  when  comparisons  are  between  cases  on  the 
same  figure.  The  choice  of  Kj/a  as  a magnification  factor  more  clearly 
shows  the  effect  of  crack  proximity  to  the  surfaces.  As  cracks  become 
larger  in  both  figures,  magnification  factors  at  0=0  are  increased  no 
more  than  15%.  In  Figure  4.5  for  cracks  with  a/c=0.5,  a/t=0.125  and 
0,1875,  the  crack  shape  effect  dominates  resulting  in  a maximum  stress 
intensity  at  about  50  degrees  of  parametric  angle  from  the  minor  axis. 
However,  when  a/t  reaches  0.25,  the  surface  effect  takes  over,  shifting 
the  point  of  maximum  stress  Intensity  to  the  hole  intersection.  The 
results  shown  in  Figure  4.6  are  similar  with  respect  to  the  strong  sur- 
face effect  but  due  to  the  larger  a/c  ratios,  the  shape  effect  is  of 
lesser  importance. 

Figure  4.7  shows  magnification  factor  versus  parametric  angle  for 
the  same  crack  with  hole  diameter  to  plate  thickness  ratios  of  1.0  and 
0.5.  As  with  similar  single  and  double  corner  crack  results,  the  re- 
duction in  d/t  involved  only  the  hole  diameter  while  all  other  plate 
and  crack  dimensions  remained  unchanged.  The  hole  stress  concentration 
effect  again  dominates  the  results  by  lowering  magnification  factors 
everywhere  since  the  crack  for  the  d/t=0.5  case  extends  into  regions  of 
lower  stress. 

4.4  Comparisons  with  Estimates  and  Solutions  of  Other  Authors 

Kobayashi  (7)  considered  the  problem  of  an  embedded  crack  emanating 
from  an  open  fastener  hole  in  a very  thick  plate.  He  used  a super- 
position method  based  upon  removing  residual  crack  surface  tractions 
predicted  by  the  stress  distribution  In  an  uncracked  plate.  The  stress 
Intensity  factor  for  a through-crack  Is  calculated  using 
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(4.4.1) 


where  Oyy(x.O)  is  the  plane  strain  stress  distribution  in  an  uncracked 
plate  with  a hole,  a"  (x,0)  is  the  mirror  image  of  o (x,0)  and  a is 

» «/  J J 

the  crack  length,  is  approximated  by  the  linear  distribution 


where  r is  the  hole  radius  and  n is  the  slope  of  the  straight  line 
distribution  which  minimizes  differences  from  the  true  distribution. 
Correction  is  made  to  the  through-crack  stress  intensity  to  produce 
corresponding  values  for  semi -elliptical  embedded  cracks. 

Figure  4.8  shows  a comparison  of  magnification  factor  versus 
parametric  angle  from  the  hole  intersection  for  an  embedded  crack  with 
a/c=0.4  as  considered  by  Kobayashi  and  a/c=0.5  from  the  finite  element- 
alternating  method.  Since  Kobayashi 's  estimate  is  for  a very  thick 
plate,  no  surface  correction  is  made.  Also,  the  plane  strain  stress 
distribution  obviously  does  not  account  for  through-thickness  stress 
variations  and  Kobayashi  admits  that  the  hole  surface  effect  could 
amount  to  a 5%  change  in  his  results. 

The  comparison  in  Figure  4.8  shows  the  magnification  factor  from 
the  finite  element -alternating  method  10?  higher  than  Kobayashi ‘s  re- 
sult at  the  hole  intersection.  At  least  5%  of  this  is  due  to  the  sur- 
face effect  unaccounted  for  in  Kobayashi 's  result.  Another  2%  can  be 
attributed  to  differences  in  describing  the  stress  distribution  over 
the  crack.  The  remaining  difference  is  accounted  for  by  the  shape  and 
hole  surface  effects,  both  of  which  have  the  tendency  to  elevate  stress 
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intensity  factors  at  the  hole  intersection  in  the  finite  element- 
alternating  method  result. 

At  0=0,  the  finite  element-alternating  method  result  is  S%  below 
that  of  Kobayashi.  The  major  portion  of  this  difference  is  from  the 
representation  of  the  stress  field  over  the  crack.  Kobayashi 's  linear 
simplification  overestimates  the  three  dimensional  Sternberg  and 
Sadowsky  (23)  solution  used  in  the  present  method  by  about  7%. 

Figure  4.9  is  a parametric  plot  of  magnification  factor  at  the 
hole  intersection  versus  2c/t  for  a series  of  results  from  the  finite 
element-alternating  method  solution  for  a/r=0.25  and  Kobayashi 's  esti- 
mate for  a/r=0.20.  The  lack  of  agreement  at  low  2c/t  is  due  to  the 
effects  mentioned  earlier.  Of  interest  is  the  large  difference,  due 
essentially  to  the  surface  effect,  as  2c/t  approaches  1.0  and  the 
nearly  constant  magnification  factor  predicted  by  the  finite  element- 
alternating  method  as  2c/t  increases  from  0.375  to  1.0. 

Browning  (3)  applied  the  finite  element-alternating  method  to  find 
magnification  factors  for  embedded  circular  cracks  emanating  from 
fastener  holes  in  finite-thickness  plates.  While  the  current  version 
of  this  method  is  not  suited  to  solution  of  the  circular  crack  problem 
as  mentioned  in  Chapter  2,  parametric  plots  show  agreement  with 
Browning's  results.  Figure  4.10  is  a plot  of  magnification  factor 
versus  2c/t  for  a/r=0.5.  As  2c/t  decreases,  the  current  study  predicts 
a reduction  in  magnification  factor  at  the  hole  intersection  due  to  re- 
duced surface  effects.  As  the  crack  approaches  a circular  shape,  how- 
ever, the  surface  effect  becomes  small  and  the  shape  effect  dominates, 
increasing  the  magnification  factor  at  this  location.  At  o=0,  the 
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reduction  in  magnification  factor  is  seen  throughout  the  range  of 
2c/t  as  expected  for  cracks  over  which  the  average  stress  has 
decreased. 


CHAPTER  5 


CRACKS  EMANATING  FROM  LOADED  FASTENER  HOLES 

5.1  Problem  Statement 

Previous  chapters  presented  stress  intensity  factor  results  for 
three  classes  of  crack  location  relative  to  an  open  hole  in  a finite- 
thickness plate.  The  plate  was  loaded  in  uniaxial  tension  perpendicular 
to  the  plane  of  the  crack.  In  the  problems  to  follow,  the  hole  is 
filled  with  a pin  or  fastener  to  which  all  or  a portion  of  the  remote 
load  is  transferred  (Figure  1.5b). 

5.2  Theoretical  Approach 

The  finite  element-alternating  method  as  applied  to  the  open 
fastener  hole  problem  assumed  symmetry  about  the  plane  z=0,  thus  only 
the  upper  portion  of  the  full  plate  was  modeled  (Figure  2.4). 
Necessarily,  the  upper  portion  of  the  hole  was  treated  as  a free  sur- 
face. If  it  can  be  shown  that  these  two  requirements  can  be  approxi- 
mately met  for  the  loaded  fastener  hole  case,  the  solution  method  need 
not  be  changed  and  the  finite  element  stiffness  matrices  from  the  open 
hole  cases  can  be  used  at  a great  savings  in  computer  time.  Crack  plane 
stress  distributions  for  the  open  hole  case  (Step  a.  Figure  2.2)  can 
then  be  replaced  by  distributions  for  partial  or  total  transfer  of  load 
to  the  filled  fastener  hole. 

5.3  Problem  Modeling  and  Assumptions 

To  approximately  meet  the  requirements  of  the  preceding  section 
and  also  to  form  a realistic  problem  the  following  assumptions  are  made: 
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1.  Fastener  Stiffness  and  Loading  Conditions.  The  fastener  is 
taken  to  be  rigid  and  is  loaded  in  double  shear  (Figure  1.5b). 

A rigid  fastener  cannot  bend  and  double  shear  loading  requires 
the  pin  to  remain  perpendicular  to  the  plate  surfaces.  Thus, 
bearing  stresses  on  the  hole  boundary  are  constant  through  the 
plate  thickness.  Resulting  normal  stresses  on  the  plane  per- 
pendicular to  the  loading  direction  at  the  minimum  section  are 
also  assumed  constant  through  the  thickness.  These  assumptions 
are  thought  to  be  realistic  since  in  practice  any  through- 
thickness variations  are  small  for  double  shear  loading. 
Fastener  rigidity  is  reflected  in  fastener  to  plate  stiffness 
ratios.  For  steel  fasteners  in  aluminum  plates  this  ratio  is 

3 to  1 while  for  the  same  fastener  in  Plexiglas,  it  is  7 to  1. 
For  the  7 to  1 ratio,  Harris,  Ojalvo,  and  Hooson  (33)  suggest 
a through-thickness  correction  of  6%  to  the  bearing  stress 
exerted  on  the  hole  by  a rigid  fastener  in  double  shear.  The 
resulting  through-thickness  variation  in  normal  stresses  on  the 
crack  plane  due  to  pin  bending  surely  is  much  smaller. 

2.  Hole  Boundary  Conditions.  The  fastener  is  assumed  to  be  smooth 
and  thus  transmits  only  radial  loads  to  the  hole.  Shear  at  the 
interface  is  assumed  to  have  little  effect  on  the  crack  plane 
normal  stress  field.  Support  for  these  assumptions  can  be 
found  in  the  results  of  Frocht  and  Hill  (30)  who  noted  only  a 
4%  decrease  in  hole  stress  concentration  factors  when  loaded 
steel  and  aluminum  fasteners  in  aluminum  plates  were  lubricated. 
Additionally,  shear  on  the  interface  does  not  affect  the  normal 
stress  distribution  at  other  locations  on  the  crack  plane  to 
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any  great  degree.  Finite  element  results  of  this  chapter  for 
loaded  smooth  pins  are  compared  to  the  analytic  solution  of 
Theocaris  (28)  which  included  shear  between  the  fastener  and 
the  hole  (Figure  5.5).  Differences  in  the  stresses  on  the 
crack  plane  in  the  vicinity  of  the  hole  predicted  by  the 
present  method  and  the  Theocaris  solution  are  small. 

3.  Fastener-to-Hole  Contact.  The  fastener  load  is  assumed  to  be 
transferred  over  the  lower  half  of  the  hole,  0£6<^tt  (Figure 
1.5b).  Harris,  et  al . , (33)  showed  that  this  contact  condition 
realistically  simulates  fastener  clearances  from  a net  fit  (no 
clearance  or  interference)  at  low  loads  to  a neat  fit  (0.4% 
clearance)  at  high  load  levels.  Greater  fastener  clearance 

or  interference  shifts  the  location  and  magnitude  of  the  maxi- 
mum tensile  stress  yet  Frocht  and  Hill  (30)  noted  for  neat  fit 
fasteners  that  this  maximum  always  occurred  on  the  hole  edge 
at  the  minimum  section. 

4,  Stress  Intensity  Factor.  Mode-one  stress  intensity  factors 
are  assumed  to  dominate  over  mode-two  effects.  Fastener  load- 
ing produces  non-zero  shearing  stresses  along  the  line  B=0 
(Figure  1.5b).  For  the  same  fastener  loading  considered  here 
but  for  double  through-cracks  (Figure  1.2),  Cartwright  and 
Ratcliffe  (47)  reasoned  that  the  mode-two  stress  intensity 
factor  was  no  more  than  1/10  the  magnitude  of  the  mode-one 
stress  intensity  factor  for  pin  loading.  Sih  (44)  studied  the 
relative  importance  of  mode-one  and  mode-two  stress  intensity 
factors  for  a through-crack  oriented  at  different  angles  to 
the  axis  of  remote  tensile  loading.  For  cracks  nearly 
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perpendicular  to  the  loading  axis,  the  mode-two  influence  on 
the  stress  intensity  factor  is  perhaps  at  the  most  2%  for  the 
ratio  suggested  by  Cartwright  and  Ratcliffe.  Therefore,  al- 
though mode-two  effects  are  present,  they  are  small  and  will 
be  neglected  in  this  work. 

5,  Displacement  Symmetry.  Displacements  in  the  x direction  along 
the  line  6=0  are  assumed  to  be  zero  (Figure  1.5b).  This  is 
obviously  not  precisely  correct.  The  asymmetry  in  loading  and 
asymptotic  behavior  of  the  normal  stresses  about  the  y-axis 
near  the  hole  cause  a displacement  and  rotation  of  this  line. 
Results  of  the  two-dimensional  finite  element  study  of  this 
chapter  show  these  displacements  and  rotations  small  enough 
ir,  all  load  transfer  cases  to  neglect. 

5.4  Calculation  of  Crack  Plane  Stress  Distributions 

The  assumptions  of  the  previous  section  are  realistic  and  at  the 
same  time  allow  the  finite  element  stiffness  matrices  formed  and  de- 
composed for  the  open  hole  cases  to  be  used  without  change  to  find 
stress  intensity  factors  for  loaded  hole  problems.  It  is  necessary, 
however,  to  describe  the  crack  plane  normal  stress  distribution  due  to 
the  loaded  fastener.  This  information  will  be  substituted  for  the  open 
hole  distribution  due  to  Sternberg  and  Sadowsky  (23)  which  is  shown  as 
Step  a in  the  finite  element-alternating  method  diagram  of  Figure  2.2. 

Crack  plane  normal  stress  distributions  for  the  loaded  fastener 
hole  cases  were  generated  using  established  two-dimensional  finite 
element  procedures  for  linear  elasticity.  This  was  necessary  for  the 
following  reasons.  Current  literature  contains  distributions  for 
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specialized  cases,  none  of  which  match  the  hole  boundary  conditions, 
plate  dimensions,  and  load  transfer  ratios  of  the  problems  considered 
here.  In  the  loaded  hole  problem  as  plate  width  changes  under  condi- 
tions of  constant  hole  diameter  and  remote  stresses,  the  stress  distri- 
bution on  the  crack  plane  varies  widely.  Thus,  it  is  imperative  that 
plate  hole  diameter  to  width  ratios  in  the  present  work  be  carefully 
matched  when  generating  crack  plane  stress  distributions  for  the  cases 
of  load  transfer  to  a fastener. 

The  finite  element  mesh  for  generation  of  crack  plane  stresses  is 
shown  in  Figure  5.1  for  d/t=1.0  where  the  hole  diameter  to  plate  width 
ratio,  X,  is  0.0833.  For  d/t=0.5  and  x=0. 04167,  the  plate  width  in 
Figure  5.1  is  simply  doubled.  The  96  elements  are  of  the  eight  node 
isoparametric  type.  Since  precise  values  for  stress  concentration 
factors  at  the  hole  were  not  an  issue,  extremely  fine  divisions  in  this 
region  were  unnecessary.  Rather,  sufficient  divisions  across  the  mini- 
mum section  to  accurately  approximate  the  gradient  was  deemed  more 
important. 

As  a test  of  accuracy,  the  plate  was  loaded  in  uniaxial  tension 
with  the  hole  unloaded.  A mesh  of  72  elements  having  six  divisions 
across  the  minimum  section  rather  than  the  eight  shown  in  Figure  5.1 
was  formed  for  a hole  diameter  to  plate  width  ratio  of  0.2.  Figure  5.2 
shows  tangential  (normal)  stress  across  the  minimum  section  versus 
distance  from  the  hole  edge.  Results  shown  on  Figure  5.2  compare  well 
with  those  from  the  analytic  solution  reported  by  Savin  (48)  for  the 
same  problem. 

For  the  loaded  fastener  hole  case,  equations  1.3.1  and  1.3.2 
describe  the  surface  tractions  applied  to  the  hole  boundary.  These 
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were  integrated  using  Gaussian  quadrature  to  produce  nodal  forces. 
Loading  conditions  were  scaled  to  produce  the  desired  load  transfer 
ratios.  Crack  plane  normal  stress  distributions  for  both  values  of 
hole  diameter  to  plate  width  and  three  load  transfer  cases  are  shown 
in  Figures  5.3  and  5,4.  Results  were  calculated  for  a plate  modulus 
of  elasticity  of  30x10^  psi  and  a Poisson's  ratio  of  0.25.  For  in- 
clusion in  the  finite  element-alternating  method,  these  distributions 
were  expressed  as  polynomials, 

5 

P(x)  = A + I Ax"  (5.4.1) 

° n=l  " 

over  a range  of  x measured  from  the  hole  edge  which  encompassed  the 
largest  crack  to  be  considered.  A least  square  fit  procedure  produced 
the  necessary  coefficients  and  good  accuracy  with  a fifth-order  poly- 
nomial. The  coefficients  and  limits  of  the  polynomial  exp.essions  are 
listed  in  Table  5.1  for  all  three  load  transfer  cases  and  the  two  hole 
diameter  to  plate  width  ratios. 

To  illustrate  the  convergent  behavior  of  this  method,  the  case  of 
all  the  remote  load  transferred  to  the  fastener  was  solved  with  meshes 
of  32,  72,  and  96  elements.  The  increase  from  32  to  72  elements 
changed  the  crack  plane  normal  stress  a maximum  of  2%  except  at  the  far 
edge  of  the  plate  where  the  difference  was  10%,  A further  increase  to 
96  elements  altered  72  element  results  less  than  1%  everywhere.  Exact 
division  of  the  element  arrangement  was  not  possible  here,  thus  this  is 
not  a strict  convergence  study. 

Figure  5.5  shows  results  due  to  the  present  work  for  72  finite 
elements  compared  to  the  analytic  solution  reported  by  Theocaris  (28) 
where  a=0.2  and  Poisson's  ratio  is  0.365  in  both  methods.  The  results 
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Table  5.1  Polynomial  Coefficients  Describing  Crack  Plane 
Pressure  for  Load  Transfer  to  a Fastener 


d/t 

d/w 

Polynomial 

Coefficients 

Q/P 

1 .0 

0.5 

0.1 

0.5 

0.04167 

^0 

18.17 

10.45 

4.272 

Ai 

-95.29 

-52.44 

-18.31 

Aa 

202.10 

111.2 

39.11 

A3 

-203.90 

-112.0 

-39.55 

A, 

97.64 

53.55 

18.95 

As 

-17.81 

-9.76 

-3.46 

(0  if  1 1.80) 

1.0 

0.0833 

Ao 

10.88 

6.96 

3.826 

Ai 

-46.43 

-28.17 

-13.56 

A2 

93.95 

57.25 

27.89 

A3 

-95.15 

-58.07 

-28.40 

A4 

46.70 

28.52 

13.97 

As 

-8.816 

-5.38 

-2.638 

(0  if  < 1.72) 

P(x)  = Aq  + 
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of  Theocaris  are  expected  to  be  somewhat  different  due  to  boundary 
conditions  at  the  hole  which  permit  shear  in  addition  to  radial  loading 
yet  great  differences  are  not  noted  near  the  hole  where  the  crack  is  to 
be  located. 

Cox  and  Brown  (32)  reported  experimental  photoelastic  stress  con- 
centration factors  at  the  minimum  section  found  by  Jessop,  Snell  and 
Hollister  for  neat  fit  loaded  fasteners.  For  the  cases  of  x=0. 04167 
and  0.0833,  the  finite  element  solution  agrees  to  within  7%' of  these 
experimental  results. 

5.5  Theoretical  Stress  Intensity  Factor  Results  for  Cracks  Emanating 

from  Loaded  Fastener  Holes 

Figures  5.6  through  5.13  present  finite  element-alternating  method 
stress  intensity  factor  results  for  cracks  emanating  from  loaded 
fastener  holes  in  finite-thickness  plates.  Results  are  shown  for  numer- 
ous crack  locations  with  respect  to  the  fastener  hole,  many  crack 
shapes,  and  two  hole  diameter  to  plate  thickness  ratios.  The  magnifi- 
cation factor  and  parametric  angle  are  as  defined  for  the  open  hole 
problems.  With  the  exception  of  Figure  5.6,  each  figure  shows  results 
for  open  holes  (Q/P=0.0)  and  100%  transfer  of  load  to  the  fastener 
(Q/P=1.0)  only.  This  is  to  facilitate  linear  superposition  to  find 
stress  intensity  factors  for  cases  of  partial  load  transfer  to  the 
fastener.  This  linear  superposition  technique  will  be  discussed  in  the 
next  section. 

All  results  for  Q/P=1,0  show  the  expected  large  increases  in 
magnification  factor,  especially  near  the  hole,  due  to  the  elevated 
stresses.  This  same  steep  stress  gradient  tends  to  magnify  subtle 
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differences  in  stress  intensity  factor  along  the  crack  periphery  seen 
in  the  open  hole  results  for  the  same  crack. 

Figures  5.12  and  5.13  show  magnification  factors  for  the  same  em- 
bedded crack  and  two  hole  diameter  to  plate  thickness  ratios.  The  open 
hole  case  was  discussed  in  Chapter  4 where  the  decrease  in  magnifica- 
tion factor  with  decreasing  d/t  was  noted.  When  this  same  situation 
occurs  in  the  loaded  hole  case,  the  asymptotic  behavior  of  the  normal 
stress  field  near  the  hole  causes  the  average  stress  in  the  region  of 
the  crack  to  increase  for  the  d/t=0.5  situation  with  respect  to  the 
d/t=1.0  case.  Thus  magnification  factors  are  now  increased  with  de- 
creasing d/t. 

Figures  5.8  and  5.9  show  results  for  single  and  double  cracks  of 
the  same  shape.  The  comments  of  Chapter  3 regarding  magnification 
factor  differences  for  single  and  double  cracks  emanating  from  an  open 
hole  are  also  true  for  the  case  of  100%  load  transfer  to  a fastener 
filled  hole.  Again,  the  maximum  difference  of  about  8%  occurs  in  the 
vicinity  of  a parametric  angle  of  30  degrees  with  smaller  differences 
occurring  at  the  hole  and  front  surface  intersections. 

5 . 6 Stress  Intensity  Factors  for  Partial  Load  Transfer  to  the  Fastener 

Stress  intensity  factors  for  partial  transfer  loads,  0.0<Q/P<1.0, 
(Figure  1.5b)  can  be  found  by  linear  superposition  of  results  for  the 
open  hole  (Q/P=0.0)  and  100%  load  transfer  (Q/P=1.0)  cases.  As  an 
example,  suppose  stress  intensity  factors  are  to  be  found  for  a single 
quarter-elliptical  crack  emanating  from  a fastener  hole  with  a/c*0.5, 
a/t=0.25,  d/t=1.0,  and  a 10%  transfer  of  load  to  the  fastener  (Q/P=0.1). 
Figure  5.14  shows  how  linear  superposition  of  the  open  hole  and  100% 
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load  transfer  results  forms  this  partial  load  transfer  and  Figure  5.6 
shows  the  magnification  factors  for  the  Q/P=0.0  and  1.0  cases.  From 
Figure  5.6,  a stress  intensity  factor  is  calculated  for  these  two  cases 
using  magnification  factors  at  some  location  along  the  crack  periphery 
and  the  adjusted  remote  stress  suggested  in  Figure  5.14.  Addition  of 
these  stress  intensity  factors  produces  the  desired  partial  load  trans- 
fer stress  intensity  factor  at  the  same  location  on  the  crack  front. 

Figure  5.6  includes  results  for  Q/P=0.1  and  0.5  which  are  not 
linear  superposition  results  but  rather  are  from  direct  calculations 
using  the  finite  element-alternating  method  and  the  appropriate  stress 
distributions  for  these  load  transfer  cases.  This  provides  a check  on 
the  accuracy  of  partial  load  transfer  stress  intensity  factors  found 
with  linear  superposition.  For  Q/P=0.5,  linear  superposition  results 
agree  to  within  2%  of  results  from  direct  finite  element-alternating 
method  calculations  and  within  4%  for  Q/P=0.1.  Thus,  Figures  5.7 
through  5.13  present  only  open  hole  and  100%  load  transfer  results 
since  linear  superposition  is  capable  of  producing  accurate  stress  in- 
tensity factors  for  any  intermediate  partial  load  transfer  to  the 
fastener.  ~ 

The  fact  that  linear  superposition  is  accurate  raises  another 
interesting  point.  Solution  2 of  the  finite  element-alternating  method 
requires  that  pressure  over  the  crack  be  expressed  as  a ten-term  poly- 
nomial (Equation  2.3.1).  For  cases  of  load  transfer  to  the  fastener, 
this  Is  accomplished  by  a least-square-fit  to  the  crack  pressure  pre- 
dicted by  the  two-dimensional  finite  element  results  of  this  chapter. 

As  the  amount  of  load  transferred  to  the  fastener  Increases,  the  normal 
stress  gradient  at  the  hole  becomes  steeper  (Figure  5.3)  and,  one  would 
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think,  subject  to  greater  errors  when  approximated  by  the  polynomial. 

On  the  contrary,  the  fact  that  linear  superposition  agrees  closely  with 
direct  calculations  for  partial  load  transfer  to  the  fastener  indicates 
that  errors  in  the  load  transfer  results  are  no  worse  than  for  the  open 
hole  case  where  the  mean  deviation  in  the  crack  pressure  averages  about 
8%  of  the  maximum  value. 

It  should  be  noted  that  this  superposition  applies  only  to  the 
hole  diameter  to  plate  width  ratios  used  in  this  work.  Theoretical 
results  were  calculated  based  upon  crack  plane  stress  distributions  for 
two  specific  plate  sizes.  The  100%  load  transfer  case  is  very  sensitive 
to  the  hole  diameter  to  plate  width  ratio  with  crack  plane  normal 
stresses  varying  greatly  with  the  changing  plate  dimensions.  While  it 
is  possible  to  use  linear  superposition  to  extend  results  to  other 
plate  widths,  it  should  be  recognized  that  errors  can  be  induced  and  in 
no  case  should  these  results  be  extended  to  hole  diameter  to  plate 
width  ratios  greater  than  0.2. 

5.7  Comparison  with  Results  of  Other  Authors 

Stress  intensity  factor  results  for  cracks  emanating  from  loaded 
fastener  holes  to  compare  with  the  present  work  are  rather  limited. 
Cartwright  and  Ratcllffe  (47)  used  experimental  compliance  methods  to 
produce  stress  Intensity  factors  for  double  through-cracks  at  a loaded 
fastener  hole  (Figure  1.2),  but  hole  diameter  to  plate  width  ratios 
are  so  different  as  to  prevent  comparison  with  the  present  work. 

Shah  (12)  applied  his  Green's  function  approach  to  estimate  stress 
intensity  factors  for  the  loaded  hole  case.  He  presents  two-dimensional 
stress  intensity  factors  for  single  (Figure  1.1)  and  double  (Figure  1.2) 
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through-cracks  and  three-dimensional  results  for  circular  cracks 
emanating  from  loaded  fastener  holes.  For  a comparison  by  the  finite 
element-alternating  method,  the  two-dimensional  through-cracks  were 
approximated  by  slender  part-elliptical  cracks  as  described  in  Chapters 
2 and  3 and  stress  intensity  factors  for  loaded  hole  cases  where 
L/R=1.5  were  calculated.  The  finite  element-alternating  method  pre- 
dicted the  stress  intensity  factor  for  the  single  through-crack  to  be 
17%  below  Shah's  result  and  for  double  through-cracks  to  be  10%  low. 

The  differing  plate  dimensions  and  resulting  expressions  for  crack 
plane  stresses  account  for  the  majority  of  these  differences.  Figure 
5.11  shows  a typical  three  dimensional  Shah  estimate  and  the  corres- 
ponding finite  element-alternating  method  result.  Shah  neglects  crack 
shape  effects  and  presents  results  for  circular  cracks  only.  This 
accounts  for  the  difference  in  magnification  factor  near  the  front  sur- 
face intersection.  Differences  near  the  hole  are  probably  due  to  the 
differing  expressions  for  the  gradient  near  the  hole  in  the  two  methods. 

Grandt  (14)  also  presents  stress  intensity  factors  for  the  single 
through-crack  at  a loaded  fastener  hole  (Figure  1.1).  A stress  dis- 
tribution due  to  Harris,  et  a1.,  (33)  for  a plane  inclined  at  6=8“ 
(Fioure  1.5b)  is  sufficiently  different  from  the  results  of  the  present 
work  that  the  stress  intensity  factor  calculated  by  the  finite  element- 
alternating  method  is  now  10%  higher  than  Grandt's  result  for  L/R=1.5. 

Comparisons  of  the  present  theory  with  a large  number  of  static 
fracture  tests  performed  by  the  author  are  presented  in  Chapter  6. 
Agreement  is  more  consistent  over  a wider  variety  of  three-dimensional 
crack  shapes  and  locations  relative  to  the  fastener  hole  than  the 
comparisons  shown  above. 


CHAPTER  6 


STATIC  FRACTURE  TESTS 


6.1  Introduction 

A series  of  static  fracture  tests  on  polymethylmethacrylate 
(PMMA)  specimens  was  performed  to  verify  the  theoretical  solution 
method  presented  in  this  work.  These  tests  included  numerous  crack 
shapes  and  locations  relative  to  bc•(■^  pen  and  loaded  fastener  holes 
in  plates  with  two  hole  diameter  to  plate  thickness  ratios  (Figures 

2.1  and  1.5).  Comparisons  with  experiments  are  made  at  locations  on 
the  crack  periphery  where  the  stress  intensity  factors  are  maximum  as 
determined  by  the  theory. 

6.2  Modeling  of  the  Fracture  Process  with  PMMA 

Many  researchers  have  used  polymers  such  as  epoxy,  polycarbonates 
and  PMMA  to  model  the  fracture  process  (21,  46,  49,  50,  51).  Trans- 
parency of  these  materials  allows  viewing  of  the  crack  front  which  is 
hidden  in  engineering  metals.  Additionally,  the  fatigue  properties  and 
fractured  surface  characteristics  of  polymers  are  similar  to  those  of 
metals. 

PMMA  is  ideally  suited  to  the  large  number  of  tests  performed  in 
the  present  experiment  due  to  its  commercial  availability  in  bulk  under 
the  trade  name  Plexiglas.  PMMA  is  also  relatively  isotropic  and 
brittle.  Under  normal  testing  conditions,  it  has  a tensile  strength 
of  about  7 X 10^  psi,  a modulus  of  elasticity  of  4.5  x 10^  psi  and  an 
elongation  at  fracture  of  about  5%. 
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6.3  Specimen  Preparation 

Specimen  blanks  were  cut  from  two  4x8  foot  sheets  of  commercial 
Plexiglas  having  a nominal  thickness  of  1/2  inch.  Each  blank  was 
oriented  in  the  same  direction  with  respect  to  the  sheet  from  which  it 
was  cut.  Specimens  were  approximately  7 inches  long  and  3 inches  wide 
where  the  length  was  limited  by  the  maximum  grip  spacing  of  the  testing 
machine.  The  width  was  chosen  to  approximate  that  of  an  infinite  plate 
for  the  open  hole  case  while  at  the  same  time  providing  economy  in  the 
number  of  specimens  fashioned  from  a Plexiglas  sheet.  Additionally, 
ell  plate  and  crack  dimensions  were  far  greater  than  the  parameter  B, 
where 

B = 2.5  . (6.2.1) 

Jones  and  Brown  (52)  suggest  that  this  should  be  the  case  for  a valid 
plane  strain  fracture  toughness  test  where  plasticity  effects  must  be 
small  and  localized. 

To  grow  single  corner  cracks  with  a/c<1.0  (Figure  2.1a),  the 
specimen  was  placed  in  bending  using  the  cracking  jig  shown  in  Figure 
6.1.  A razor  blade  was  touched  to  the  tension  side  of  the  blank  pro- 
ducing a semi-elliptical  surface  crack  whose  center  is  offset  the 
length  of  a hole  radius  from  the  specimen  centerline.  The  fastener 
hole  and  remote  loading  holes  were  then  carefully  drilled  and  reamed. 

Preparation  of  double  crack  specimens  with  a/c<1.0  (Figure  2.1b) 
was  similar  with  one  important  exception.  Here,  the  dual  cracks  were 
grown  in  the  same  cracking  Jig  using  a specially  fashioned  razor  blade 
which  simultaneously  produced  separate  surface  cracks  on  opposite 
sides  of  the  hole  location.  As  before,  the  fastener  hole  was  then 
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drilled  and  reamed  such  that  the  major  axis  of  each  surface  crack  was 
bisected  by  the  hole. 

A new  device  was  designed  and  procedures  developed  to  grow  single 
and  double  cracks  with  a/c>1.0  (Figures  2.1a  and  2.1b)  and  embedded 
cracks  (Figure  2.1c)  emanating  from  fastener  holes.  First,  the  fastener 
holes  and  appropriate  loading  holes  were  drilled  and  reamed.  The  speci- 
men was  then  loaded  at  remote  holes  offset  from  the  specimen  centerline 
to  produce  combined  axial  loading  and  bending.  The  crack  was  introduced 
on  the  side  of  the  fastener  hole  where  the  tensile  stresses  due  to  the 
axial  loading  were  reduced  by  the  compression  due  to  bending.  This 
allowed  for  better  control  and  prevented  unstable  growth  during  place- 
ment of  the  crack. 

The  cracking  jig  is  shown  in  place  on  a specimen  in  Figure  6.2c. 

A small  razor  blade  is  afixed  to  a rod  which  is  aligned  parallel  to 
the  hole  centerline.  With  load  on  the  specimen,  the  jig  is  tightened 
which  presses  the  blade  against  the  edge  of  the  hole.  Centering  the 
jig  at  the  center  of  the  plate  through-the-thickness  produces  an  em- 
bedded crack  while  offset  of  the  jig  toward  a plate  surface  produces  a 
single  corner  crack  at  the  hole.  Double  corner  cracks  were  grown  in  a 
similar  manner.  First  one  crack  was  placed  as  explained  above  for  the 
single  deep  corner  crack.  The  load  was  then  reduced  and  the  cracking 
jig  rotated  to  the  opposite  side  of  the  plate  where  the  second  crack 
was  grown. 

In  all  cases,  crack  sizes  and  shapes  could  be  controlled  with 
razor  blade  length  and  loading  conditions.  Exact  duplication  of  crack 
geometries  was,  however,  not  possible  yet  the  growth  of  cracks  with 
similar  geometries  was  deemed  important  due  to  the  Inherent  scatter 
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in  fracture  tests  of  this  type.  Thus,  comparisons  with  the  present 
theory  are  always  made  with  more  than  one  experimental  result. 

6.4  Slow  Growth  in  PMMA 

Fracture  in  PMMA  frequently  includes  slow  growth  preceding  rapid 
crack  extension.  The  region  of  slow  growth  on  the  fractured  surface 
is  rippled  while  the  region  of  rapid  crack  extension  is  identified  by 
a mirrored  surface.  This  is  clearly  seen  on  all  the  fractured  sur- 
faces pictured  in  Figure  6.3.  In  tests  performed  here,  slow  growth 
began  at  specific  points  on  the  crack  front  which  could  be  identified 
visually.  The  location  of  initial  slow  growth  depended  upon  the  crack 
shape  and  orientation  with  respect  to  the  hole  and  can  be  assumed  to 
coincide  with  the  location  of  maximum  stress  intensity.  The  benefit 
here  is  that  experimental  and  theoretical  comparisons  can  be  done  at 
a specific  location  on  the  crack  periphery  rather  than  the  gross  com- 
parison necessary  if  the  experiment  does  not  show  where  the  maximum 
stress  intensity  occurs. 

A series  of  load  versus  crack  opening  displacement  curves  was 
obtained  for  surface  flawed  PMMA  specimens  with  and  without  fastener 
holes.  This  was  done  in  order  to  assess  the  relative  importance  of 
the  slow  growth  and  the  accuracy  of  the  visual  procedure  in  predicting 
its  onset.  A small  linear  differential  transformer  was  mounted  over 
the  cracks  and  connected  to  an  x-y  plotter.  A typical  curve  is  shown 
In  Figure  6.4  for  double  surface  cracks  emanating  from  a fastener  hole. 
In  all  cases,  the  onset  of  slow  growth  could  be  visually  seen  to  begin 
at  loads  no  more  than  3%  above  those  predicted  by  the  load-displacement 
records.  The  load  at  the  onset  of  slow  growth,  determined  visually. 
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was  taken  as  the  critical  load  for  calculation  of  an  experimental  mag- 
nification factor  to  be  compared  with  the  theoretical  predictions.  The 
location  of  the  initial  growth  was  assumed  to  be  the  location  of 
maximum  stress  intensity. 

A series  of  experiments  was  recorded  on  high-speed  16  mm  movie 
film.  Analysis  of  this  film  qualitatively  supports  the  observations 
noted  above.  Additionally,  when  these  films  are  used  in  conjunction 
with  inspection  of  the  fractured  surfaces  and  parametric  plots  of 
stress  intensity  factors  from  the  present  theory,  interesting  conclu- 
sions regarding  crack  growth  can  be  formed.  The  fractured  surfaces 
from  some  of  the  movie  tests  are  pictured  in  Figure  6.3  and  remarks  on 
the  slow  growth  seen  will  be  made  later  in  this  chapter. 

6.5  Specimen  Testing 

Prior  to  testing,  specimen  widths  and  thicknesses  were  recorded. 
All  specimens  were  loaded  to  failure  in  an  Instron  Testing  Machine 
at  a loading  rate  of  0.02  inches  per  minute  to  avoid  strain  rate 
effects.  For  each  specimen,  the  location  on  the  crack  periphery  at 
which  slow  growth  began  and  the  corresponding  load  were  noted.  Crack 
dimensions  were  measured  with  a dial  gage  following  failure. 

The  loading  arrangements  for  open  holes  and  100%  load  transfer  to 
a fastener  are  shown  in  Figures  6.2a  and  6.2b.  All  loading  fixtures 
were  carefully  machined  and  aligned  to  eliminate  plate  bending  effects. 
The  fastener  was  machined  to  0.4%  clearance  to  simulate  a neat  fit  and 
lubrication  was  applied  to  approximate  the  smoothness  condition  at  the 
hole-fastener  interface. 
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The  laboratory  atmosphere  was  not  controlled  specifically  during 
testing.  Minor  changes  in  temperature  and  humidity  which  did  occur 
are  deemed  not  significant. 

6.6  Data  Reduction 

A data  reduction  method  referred  to  as  the  baseline  method  is  used 
in  this  work.  This  method  has  been  applied  in  a number  of  earlier 
studies  (46,  50,  53)  and  will  be  used  here  as  follows.  For  a semi- 
elliptical surface  crack  in  a finite  thickness  plate,  the  stress 
intensity  factor  can  be  defined  as 

Kj  = Ml  M2  o/WQ  (6.6.1) 

where  Mi  and  M2  are  the  front  surface  and  back  surface  correction 
factors  respectively,  and  Q is  the  usual  crack  shape  parameter.  The 
front  surface  correction  factor.  Mi,  is  a function  of  a/2c  while  the 
back  surface  correction  factor,  M;,  is  a function  of  both  a/2c  and  a/t. 
As  a/t  approaches  zero,  the  back  surface  effect  disappears  and  M2 
approaches  a value  of  unity.  Smith  and  Alavi  (37)  solved  for  stress 
intensity  factors  around  the  periphery  of  a semi -elliptical  crack  in 
the  surface  of  a semi -infinite  solid.  For  small  a/2c  values,  the 
deepest  position  on  the  crack  front  is  most  critical  and  the  above 
authors  approximated  the  front  surface  correction  factor  at  this  loca- 
tion by 

M - 1.12  - .0476  (a/2c)  - .227  (a/2c)2,  (a/2c  i0.42)  . (6.6.2) 

••or*  r*cent  Investigations  by  Smith  and  Sorensen  (39)  and  Shah  and 

(n)  predict  front  surface  effects  which  are  about  2%  lower 
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over  the  range  of  crack  shapes  considered  here.  Nevertheless,  Equation 
6.6.2  provides  a useful  functional  relationship  for  the  purposes  of 
this  study.  Thus  for  low  values  of  a/2c  and  small  a/t,  the  front  sur- 
face correction  is  known  and  the  back  surface  correction  can  be  taken 
as  unity.  If  the  remote  failure  stress  and  crack  dimensions  can  be 
measured  using  suitable  surface  cracked  baseline  specimens,  the 
critical  value  of  Kj,  or  the  fracture  toughness,  can  be  determined  for 
the  material.  This  value  will  be  defined  as  the  apparent  fracture 
toughness , 

= Ml  a»WQ'  . (6.6.3) 

Separate  baseline  tests  were  conducted  for  each  of  the  two  PMMA 
sheets.  Baseline  specimens  were  approximately  1.5  inches  wide  versus 
3 inches  for  the  fastener  hole  specimens.  All  other  dimensions  were 
the  same  as  for  the  fastener  hole  tests.  Smith  (46)  showed  that  for 
surface  flaws,  specimen  width  effects  were  unimportant  for  a/2c  < 0.3 
if  w/2c  is  larger  than  3 where  w is  the  specimen  width.  Several  base- 
line specimens  fell  slightly  outside  this  range  yet  no  evidence  of 
width  effects  could  be  seen. 

The  results  of  the  two  baseline  studies  are  shown  in  Figures  6.5 
and  6.6.  The  values  of  AKj^,  were  taken  as  the  average  of  values  for 
all  specimens  tested  from  each  plate.  The  maximum  data  scatter  in  each 
test  was  ± 7%.  Fracture  toughness  results  obtained  here  agree  well 
with  those  of  other  experiments  (55). 

The  slow  growth  phenomena  in  PMMA  raises  questions  about  the 
validity  of  the  apparent  fracture  toughness  va*lues  obtained  from  the 
baseline  tests  performed  here.  Brown  and  Srawley  (56)  commented  to 
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great  lengths  on  the  original  concept  of  Irwin  (57)  regarding  resist- 
ance to  crack  extension  during  a fracture  toughness  test.  Irwin  pro- 
posed that  the  crack  extension  force,  G,  was  balanced  by  increasing 
resistance  to  crack  extension,  R.  A condition  of  equilibrium  between 
G and  R is  maintained  to  the  point  of  crack  instability.  Should  large 
amounts  of  crack  slow  growth  occur,  it  is  possible  that  the  instability 
point  could  be  located  incorrectly,  resulting  in  values  for  fracture 
toughness  which  are  in  error.  Crack  opening  displacement  plots  made  on 
surface  flawed  baseline  specimens  in  this  work  showed  slow  growth  to 
occur  within  no  more  than  of  the  maximum  load.  Additionally,  this 
growth  was  very  rapid  and  nearly  impossible  to  view,  being  detected 
primarily  on  the  fractured  surfaces  after  failure.  Therefore,  the 
author  feels  slow  growth  in  the  PMMA  baseline  tests  performed  here  to 
have  little  effect  on  the  resulting  fracture  toughness  values  presented. 
Also,  since  the  load  at  the  onset  of  slow  growth  and  the  load  at  which 
specimen  separation  occurs  are  so  close,  the  separation  load  is  used  to 
calculate  the  apparent  fracture  toughness  values. 

The  apparent  fracture  toughness  found  from  the  baseline  tests  is 
used  in  the  reduction  of  experimental  data  from  the  fastener  hole  tests 
as  follows.  The  magnification  factor  is  defined  as 


o/irF 

where  a is  the  remote  stress  at  the  onset  of  slow  growth  in 
hole  test  and  b is  the  crack  minor  axis  regardless  of  crack 
For  cases  of  100%  load  transfer  to  the  fastener,  the  stress 
ed  as  a bearing  stress. 


(6.6.4) 

an  open 
orientation, 
is  express- 
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where  P is  the  remote  load  at  the  onset  of  slow  growth,  d is  the  hole 
diameter,  and  t is  the  plate  thickness.  This  is  necessary  in  compar- 
isons with  the  theory  due  to  differing  theoretical  and  experimental 
plate  widths.  The  definition  of  magnification  factor  shown  above  does 
not  separate  the  various  effects  on  the  stress  intensity  factor  men- 
tioned in  Chapter  2 but  describes  a combined  effect  for  each  problem. 
This  is  entirely  suitable  for  the  purpose  of  verification  of  the 
present  theory  since  the  theoretical  magnification  factors  predicted 
in  this  work  also  do  not  distinguish  between  these  effects. 


6.7  Comparison  of  Theory  with  Experiments 

Figures  6.7  through  6.22  present  comparisons  of  experimental  and 
theoretical  magnification  factors.  Each  figure  shows  results  for  a 
particular  crack  location  with  respect  to  open  and  loaded  fastener 
holes.  Experimental  magnification  factors  as  defined  in  the  previous 
section  are  plotted  versus  the  non-dimensional  crack  depth,  a/t,  and 
the  range  of  crack  shapes  in  each  figure  is  small.  Theoretical  calcu- 
lations were  made  for  crack  shapes  and  depths  which  are  close  to  the 
average  experimental  crack  shape  and  depth  values  for  the  experiments 
plotted.  The  theoretical  magnification  factor  is  defined  as 


(6.7.1) 


for  open  hole  cases  or 


ajj/nF 


(6.7.2) 
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for  100%  transfer  of  load  to  a fastener  where  and  b are  as  defined 
for  the  experimental  magnification  factor.  Problem  dimensions  are  as 
depicted  in  Figure  2.1. 

Magnification  factor  comparisons  are  made  at  points  on  the  crack 
periphery  corresponding  to  the  location  of  onset  of  slow  growth  in  the 
experiments.  In  most  cases  this  also  corresponded  to  the  location  of 
maximum  theoretical  magnification  factor.  One  exception  is  the  series 
of  tests  for  hole  diameter  to  plate  width  ratios,  d/t,  of  0.5.  These 
specimens  exhibited  little  or  no  slow  growth  so  it  was  impossible  to 
determine  from  the  experiment  at  what  position  along  the  crack  border 
to  make  comparison  with  the  theory.  Therefore,  comparisons  are  made  at 
locations  on  the  crack  periphery  where  the  theory  predicts  the  magnifi- 
cation factor  to  be  maximum.  Another  exception  is  the  embedded  crack 
where  the  slow  growth  region  was  very  small,  but  it  covered  a range  of 
the  crack  border  encompassing  the  location  of  the  maximum  theoretical 
magnification  factor. 

Experimental  data  scatter  in  11  of  the  15  comparative  plots  is 
less  than  ± 9%.  The  remaining  four  figures  (6.9,  6.13,  6.16,  and  6.20) 
show  larger  data  scatter.  This  is  generally  the  result  of  one  data 
point  being  separated  from  an  otherwise  close  grouping  of  experimental 
results.  Scatter  in  these  four  figures  is  reduced  to  less  than  ± 10% 
if  the  separated  data  points  are  not  considered.  However,  if  all 
points  are  used,  experimental  data  scatter  ranges  from  ± 11%  to  ± 18% 

In  these  four  cases. 

Theoretical  calculations  are  for  a material  with  a Poisson's  ratio 
of  0.25  while  for  PMMA  this  value  is  about  0.35.  Smith  (46)  noted  a 
3%  increase  in  theoretical  back  surface  correction  factors  with  an 


69 


increase  in  Poisson's  ratio  from  0.25  to  0.39  in  studies  of  semi- 
elliptical surface  cracks.  Assuming  this  increase  would  also  apply  in 
the  present  work,  all  theoretical  magnification  factors  should  be 
raised  3%.  However,  it  was  noted  in  an  earlier  section  that  the  visual 
procedure  used  to  predict  loads  at  the  onset  of  slow  growth  in  PMMA  was 
about  3%  high.  If  such  a correction  was  applied  to  all  experimental 
points,  magnification  factors  would  increase  3%.  Since  the  Poisson's 
ratio  effect  and  the  correction  to  the  visual  procedure  both  predict 
adjustment  in  the  same  direction,  all  experimental  and  theoretical 
points  on  the  figures  would  simply  be  higher  by  about  3%.  The  off- 
setting nature  of  these  corrections  is  rather  unimportant  in  the  present 
comparison  and  thus  is  not  shown  on  the  figures  to  follow. 

For  purposes  of  explaining  the  differences  between  theoretical  and 
experimental  results  shown  in  Figure  6.7  through  6.22,  a percent  differ- 
ence will  be  defined  as: 

% Difference  = ' "'^Experimental  ^ ,5  7 3, 

^^Theoretical 

Figure  6.7  through  6.10  show  comparisons  of  experimental  and 
theoretical  magnification  factors  for  single  corner  cracks  adjacent  to 
both  open  and  loaded  fastener  holes  where  the  hole  diameter  to  plate 
width  ratio,  d/t,  is  1.0.  The  open  hole  comparisons  are  within  8% 
while  the  two  cases  of  100%  load  transfer  to  the  fastener  show  differ- 
ences of  13%  and  17%. 

Double  crack  comparisons  for  d/t=1.0  and  both  open  and  loaded 
holes  are  shown  in  Figures  6.11  through  6.14.  The  open  hole  comparison 
In  Figure  6.11  is  very  close  while  in  Figure  6.12  the  average  experi- 
mental value  is  about  6%  below  the  theoretical  prediction.  Figures 
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6.13  and  6.14  show  100%  load  transfer  comparisons  and  have  the  experi- 
mental results  about  16%  and  5%  lower,  respectively,  than  the  theory. 

Comparisons  of  theoretical  and  experimental  magnification  factors 
for  embedded  cracks  are  shown  in  Figures  6.15  and  6.16  for  d/t=1.0. 

The  open  hole  case  of  Figure  6.15  has  the  experimental  results  about 
9%  low  while  the  100%  load  transfer  case  comparison  in  Figure  6.16  is 
very  good. 

Figures  6.17  through  6.22  show  comparisons  for  cases  of  d/t=0.5. 

All  three  crack  locations  with  respect  to  the  hole  are  represented. 

With  the  exception  of  Figure  6.22  which  is  for  100%  load  transfer  to 
a fastener,  comparisons  shown  are  for  open  holes.  The  embedded  crack 
cases  of  Figures  6.21  and  6.22  show  very  good  agreement  with  the 
present  theory.  Single  and  double  crack  comparisons  in  Figures  6.17 
through  6.20  have  the  experimental  results  from  19%  to  31%  lower  than 
the  theory. 

Although  the  agreement  between  theoretical  and  experimental  mag- 
nification factors  in  most  cases  presented  here  is  good,  the  differences 
which  do  occur  can  perhaps  be  at  least  qualitatively  explained.  The 
existence  of  a small  yield  zone  around  the  periphery  of  a crack  has 
been  well  established.  The  size  of  this  crack-tip  yield  zone  increases 
with  increasing  local  stresses  in  the  neighborhood  of  the  crack.  All 
cracks  in  this  study  are  located  in  a region  of  high  stress  concentra- 
tion near  a hole  in  a plate.  Additionally,  comparisons  between  experi- 
ments and  theory  are  generally  made  at  points  of  Intersection  of  the 
crack  with  a plate  surface.  Crack-tip  plastic  zones  at  these  locations 
are  expected  to  be  larger  than  at  other  locations  along  the  crack 
periphery. 
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Using  experimental  failure  loads,  estimates  of  stresses  in  the 
vicinity  of  cracks  emanating  from  fastener  holes  in  the  PMMA  tests 
performed  here  can  be  made.  It  is  then  possible  to  obtain  a quali- 
tative judgment  of  the  importance  of  plastic  deformation  in  the  experi- 
ment. A plot  of  percent  deviation  between  theory  and  experiment  versus 
the  ratio  is  shown  in  Figure  6.23  for  all  cases  except  embedded 

cracks  which  will  be  discussed  later.  A strong  trend  in  this  plot  in- 
dicates deviations  from  the  theory  of  from  12%  to  31%  for  experiments 
in  which  exceeded  0.75.  However,  when  was  less  than 

0.75,  agreement  was  good  and  the  deviation  never  exceeded  8%.  It  is 
thus  not  possible  to  ignore  plasticity  effects  in  the  experiments  for 
high  ratios  of  o^g^^^ys  probably  accounts  for  much  of  the  devia- 

tion from  the  theory  for  the  d/t=0.5  and  100%  load  transfer  cases.  ^ 

The  closest  theoretical  to  experimental  agreement  for  a particular 
crack  location  with  respect  to  a fastener  hole  occurred  for  cases  of 
embedded  cracks  regardless  of  loading  conditions  or  hole  diameter  to 
plate  thickness  ratios.  Although  local  stresses  were  also  high  in 
these  tests,  magnification  factors  are  compared  at  a location  well  re- 
moved from  the  hole  surface  and  plate  lateral  surfaces.  Yield  zones 
are  most  probably  smaller  at  this  location  than  at  surface  intersections 
where  comparisons  are  made  for  most  other  cases. 

6.8  Some  Observations  on  Crack  Growth  in  the  Present  Experiments 
The  experiments  conducted  by  the  author  were  for  the  expressed 
purpose  of  verification  of  the  present  theory.  In  the  course  of  the 
experimentation,  however,  interesting  slow  growth  near  cracks  emanating 
from  fastener  holes  in  PMMA  was  noted.  For  certain  cases,  this 
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phenomena  was  recorded  on  movie  film  and  for  all  tests  the  fractured 
surfaces  were  carefully  inspected.  An  analysis  of  the  growth  behavior 
using  the  present  theory  is  of  sufficient  interest  to  be  included  here. 

Figure  6.3  shows  photographs  of  six  specimens  following  fracture. 
High  speed  16  mm  movies  were  also  taken  of  these  same  specimens  during 
loading.  The  upper  four  photographs  in  the  figure  are  open  hole  cases 
and  have  the  initial  crack  shape  outlined  for  better  definition.  The 
bottom  two  pictures  are  for  100%  transfer  of  load  to  the  fastener. 

Here,  the  initial  crack  shape  is  more  easily  seen  and  is  not  outlined. 
The  f«^actured  surfaces  on  the  six  specimens  are  representative  of  sur- 
faces seen  in  all  experiments. 

In  the  photographs  of  Figure  6.3,  the  regions  of  slow^ growth  are 
the  rippled  surfaces  while  rapid  crack  extension  is  typified  by  the 
mirrored  surfaces.  Theocaris  (59)  explained  that  the  wavy  surface  of 
the  slow  growth  region  is  caused  by  a similarly  wavy  elastic  area  out- 
side a circular  highly  strained  zone  around  the  crack  tip.  As  the 
crack  grows  through  this  rippled  elastic  region,  the  fractured  surface 
likewise  becomes  wavy. 

Figures  6.3a  and  6.3b  show  fractured  surfaces  for  single  and 
double  corner  cracks  of  low  a/c.  Slow  growth  initiates  at  the  hole 
intersection  in  both  cases  and  growth  patterns  are  similar.  Rarely  is 
any  growth  seen  at  the  front  surface  Intersection.  Figure  6.24  is  a 
plot  of  Kj/a  versus  a/c  using  the  present  theoretical  results  In  an 
attempt  to  explain  this  growth  pattern  where  slow  growth  is  depicted 
from  left  to  right  on  the  figure.  As  slow  growth  proceeds  from  regions 
near  the  hole,  stress  intensity  factors  increase  at  both  the  hole  and 
surface  intersections  but  more  so  at  the  front  surface.  The  crack 
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grows  in  such  a fashion  as  to  attempt  to  equalize  stress  intensity 
factors  all  along  the  crack  periphery. 

Figure  6.3c  shows  a photograph  of  double  cracks  which  are  longer 
in  the  hole  direction  than  along  the  front  surface.  Here  and  in  similar 
single  crack  cases,  slow  growth  begins  at  the  front  surface  inter- 
section but  soon  occurs  everywhere  along  the  crack  periphery.  Figure 
6.25  is  a parametric  plot  of  theoretical  results  for  this  observed 
growth  pattern  where  slow  growth  is  depicted  from  left  to  right  on  the 
figure.  At  step  1 on  this  figure  (a/c=3.0,  a/t=0.5)  the  stress  inten- 
sity factor  at  the  front  surface  intersection  is  largest  so  growth 
begins  there.  Between  Steps  1 and  2 (a/c=2.0,  a/t=0.5),  the  location 
of  maximum  stress  intensity  shifts  to  the  hole  intersection  and  growth 
begins  at  this  point  also.  Steps  3 and  4 (a/c=1.5,  a/t=0.75,  0.9)  show 
increases  everywhere  as  growth  continues  all  along  the  crack  border. 

The  fractured  surface  shown  in  Figure  6.3d  is  typical  for  an 
embedded  crack  emanating  from  an  open  hole  in  PMMA.  A small  slow 
growth  area  is  seen  between  the  initial  crack  and  one  plate  lateral  sur- 
face. The  remaining  fractured  surface  is  mirrored.  Figure  6.26  is  a 
parametric  plot  of  theoretical  results  approximating  this  case  where 
slow  growth  is  depicted  from  right  to  left  on  the  figure.  Growth  is 
assumed  to  initiate  at  the  location  of  maximum  stress  intensity,  in  this 
case  at  a parametric  angle  of  about  40  degrees  from  the  crack  minor 
axis.  The  crack  shape  changes  in  such  a manner  that  stress  intensities 
rapidly  increase  at  this  location  but  decrease  at  the  hole  intersection 
and  rapid  ^racture  follows. 

Figures  6.3e  and  6.3f  are  photographs  of  typical  fractured  sur- 
faces for  two  cases  in  which  100-li  of  the  remote  load  is  transferred  to 
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the  fastener-filled  hole.  Here,  either  rapid  crack  extension  or  slow 
growth  followed  by  rapid  crack  extension  begins  at  the  hole  intersec- 
tion. Crack  arrest  occurs  when  the  crack  shape  approximates  that  of  a 
through-crack  for  which  stress  intensity  factors  are  lower.  This  is 
followed  by  more  slow  growth  of  the  new  through-crack  shape  over  rather 
large  distances  across  the  plate  half-width. 
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I • CHAPTER  7 

I 'CRACK  OPENING  DISPLACEMENTS 

m 7.1  Problem  Statement 

* Unlike  the  method  presented  in  this  work,  some  stress  intensity 

I factor  solution  use  prior  knowledge  of  crack  opening  displacements. 

The  calculation  of  these  displacements,  while  not  a primary  objective 
I of  this  work  since  stress  intensity  factors  are  calculated  indepondent- 

_ ly,  is  sometimes  useful  in  extending  solutions  to  cracks  of  other 

" geometries.  Therefore,  a simple  relationship  for  opening  displacements 

I of  surface  cracks  emanating  from  fastener  holes  is  presented  in  this 

section. 

^ 7.2  Solution  Method 

■ One  of  the  solutions  used  in  the  finite  element-alternating  method 
is  that  of  an  embedded  elliptical  crack  in  an  infinite  solid.  Shah  and 

I Kobayashi  (13)  extended  the  solution  by  Green  and  Sneddon  (41)  for  this 

m problem  to  allow  arbitrary  normal  loading  on  the  crack  where  this 

loading  is  expressible  as  a ten-term  polynomial  in  x and  y.  The  details 
/ of  their  work  are  widely  known  and  will  not  be  presented  here.  However, 

a brief  description  of  the  problem  and  their  solution  method  is  neces- 
I sary  to  provide  background  for  the  present  work. 

m Navier's  displacement  equations  of  elasticity  are  identically 

satisfied  by  a harmonic  function,  ij),  if  ^ satisfies  Laplace's  equation, 

■ 7^^  ^ 0 , (7.2.1) 

I provided  shear  is  zero  on  the  plane  z = 0.  The  z displacement  com- 

f ponent  of  interest  here  is 
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w = z|^-2(l-,)||.  . 


(7.2.2) 


F 
I 


On  the  crack  plane,  z = 0,  this  equation  reduces  to 


w = - 2(l-n)  If 


(7.2.3) 


where  n is  Poisson's  ratio. 

Ellipsoidal  coordinates  x,  p,  and  v are  introduced  to  solve  the 
boundary  value  problem,  where  these  coordinates  are  roots  of  the 
equation  in  s. 


-1=0, 

a +s  b +s  s 


(7.2.4) 


and  the  ellipsoidal  and  cartesian  coordinates  are  related  by 

a2(a2-b2)x2  = (a2+x)(a2+p)(a2+v)  , (7.2.5) 

b2(b2-a2)y2  = (b2+x)(b2+p)(b2+v)  , (7.2.6) 

a2b2z2  = Apv  , (7.2.7) 

subject  to  the  standard  restriction, 

-a2£vj^-b2^p£0^X<“ 

In  the  ellipsoidal  coordinate  system,  x = 0 represents  the  region 
in  the  z = 0 plane,  interior  to  the  crack,  while  p = 0 represents  the 
region  outside  the  crack  in  this  same  plane. 

Segedin  (60)  suggested  the  use  of  a harmonic  potential  function 


y(n)  _ J w'^(s)ds 

X mr 


(7.2.8) 
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where 


w(s)  = 1 


(7.2.9) 


and 


Q(s)  = s(a2+s)(b2+s) 


(7.2.10) 


The  general  stress  function,  ij),  is  formed  as  a polynomial  in  V 
3 3 


(n) 


- I I 

i=0  j=0 


_LL 


(i+j<3) 


(7.2.11) 


or  i)>  itself  can  be  expressed  as  the  double  series, 


3 3 

4.  = I I (i+Jl3)  . (7.2.12) 

i=0  j=0 

From  equation  (7.2.3),  the  calculation  of  crack  opening  displacements 
involves  only  the  first  derivative  of  <|)  and  each  term  in  the  series 
(7.2.12)  can  be  differentiated  separately.  In  the  interest  of  brevity, 
only  will  be  presented  in  detail  to  illustrate  the  procedure. 

From  equation  (7.2.8), 


- r 7 w(s)dS 

1 /(JITT 

where  0^^  differs  from  c^by  an  integral  constant, 
rule  for  differentiating  integrals. 


(7.2.13) 

Using  the  Leibnitz 


9*^  CO 

IT  - ■ (7.2.14) 

X s/?)rsy 

The  integral  in  equation  (7.2.14)  can  be  expressed  in  terms  of  elliptic 
integrals  and  Jacobian  elliptic  functions  and  the  limiting  form  of  ■— 
can  be  found  by  setting  y = 0 and  approaching  the  crack  tip  from 
outside.  While  this  seems  straightforward,  the  work  necessary  to 
obtain  the  ten  first  derivatives  needed  for  crack  opening  displacements 
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is  excessive  and  tedious.  Nevertheless,  several  derivatives  were  cal- 
culated  in  this  fashion  resulting  in  answers  identical  to  those  obtained 
by  the  easier  approach  which  follows. 

If  a limiting  process  in  z is  applied  directly  to  the  integral  in 
equation  (7.2.14),  the  amount  of  tedious  calculation  and  bookkeeping  is 
considerably  reduced.  Approaching  the  crack  plane  from  above. 


8$ 


(7.2.15) 


where 

f(z)  = / . (7.2.16) 

X s/Q(s) 

Since  equation  (7.2.15)  is  an  indeterminate  form,  applying  L'Hospitals 
rule  yields 

3 ^ \ 

JT  ^ 

Again,  using  the  Leibnitz  rule. 


f'(z)  = - 


x/QTxT 


11 

dz 


(7.2.18) 


Calculating  (7.2.7),  substituting  Q(x)  from  (7.2.10),  and 

performing  the  necessary  algebra  leads  to 


3$ 

” _ 00 
3i FFT- 


(7.2.19) 


where  w = a^b^  - b^x^  - a^y^  , The  calculation  of  the  remaining 
nine  derivatives  proceeds  in  a similar  fashion.  The  final  crack  open- 
ing displacement  equation  is  shown  in  Appendix  B. 

To  check  the  limiting  process  used  in  obtaining  the  ten  derivatives 
along  with  the  algebra  involved,  crack  opening  displacements  were 
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calculated  for  a circular  crack  under  constant  and  linearly  varying 
normal  pressure.  The  circular  crack  was  approximated  by  a near-circular 
elliptical  crack  with  a/b  = 0.99.  Since  this  analysis  is  being  done 
algebraically,  there  is  no  concern  about  the  ratio  a/b  causing 
numerical  difficulties  such  as  those  mentioned  in  Chapter  2.  Shah  and 
Kobayashi  (22)  provide  the  necessary  stress  function  coefficients  for 
the  two  loading  cases.  Smith  (45)  solved  this  same  problem  for  loading 
expressed  in  the  form  of  a Fourier  series.  Crack  opening  displacements 
calculated  by  the  present  method  compare  with  Smith's  values  to  within 
0.50%  for  both  loading  cases.  This  small  difference  is  probably  due  to 
the  approximation  of  a circular  crack  with  a near-circular  elliptical 
crack. 

7.3  Solution  Results 

A necessary  prerequisite  to  the  calculation  of  crack  opening  dis- 
placements using  the  method  of  this  section  is  a knowledge  of  the 
stress  function  coefficients,  C..,  for  a given  problem  geometry.  For 

* V 

each  iteration  of  the  finite  element-alternating  method,  stress  in- 
tensity factors  and  stress  function  coefficients  are  calculated  based 
upon  the  freeing  of  normal  stress  on  the  crack  plane. 

A final  set  of  stress  function  coefficients  for  a particular  case 
is  formed  by  sunning  like  coefficients  for  each  iteration.  Appendix  A 
lists  the  final  stress  function  coefficients  for  all  problems  solved. 
Crack  opening  displacements  are  found  by  using  the  equation  in 
Appendix  B together  with  a set  of  stress  function  coefficients  from 
Appendix  A.  Figure  7.1  shows  a typical  result  In  the  form  of  a three- 
dimensional  plot  of  crack  opening  displacements  for  the  case  of  a 
single  corner  crack  emanating  from  an  open  hole  In  an  Infinite  plate. 
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CONCLUSION 

8.1  Discussion 

A new  method  developed  by  Browning  (3)»  Ganong  (4)  and  the  author 
was  used  to  predict  stress  Intensity  factors  for  three-dimensional 
linear  fracture  mechanics  problems  not  previously  solved.  This  tech- 
nique, called  the  finite  element-alternating  method,  produced  stress 
Intensity  factors  around  the  periphery  of  surface  cracks  emanating  from 
fastener  holes  In  finite- thickness  plates.  Problems  solved  Included 
three  crack  locations  with  respect  to  the  hole  and  numerous  crack 
shapes.  Two  hole  diameter  to  plate  thickness  ratios  were  considered 
and  results  were  presented  for  both  open  and  fastener-filled  ' to 
which  a portion  of  the  remote  load  Is  transferred. 

The  solution  method  presented  Is  accurate,  economical  and  capable 
of  modeling  complex  three-dimensional  geometries.  Compar1son'>  of 
present  results  with  estimates  of  other  authors  showed  similar  trends 
but  better  accuracy  due  to  the  accounting  for  all  effects  In  the  prob- 
lems considered.  Agreement  with  Independent  fatigue  experiments  by 
Snow  (20)  Is  good  for  the  several  crack  shapes  and  depths  considered. 
Comparisons  of  the  present  theory  with  stress  freezing  photoelastic 
experiments  by  McGowan  and  Smith  (19)  are  erratic  at  best.  The  author 
feels  that  the  results  from  these  experiments  are  In  error,  particularly 
at  locations  of  crack  Intersection  with  a hole. 

Static  fracture  tests  In  PMNA  were  conducted  to  provide  additional 
verification  of  the  theoretical  calculations  performed.  New  experimen- 
tal procedures  allowed  for  testing  of  specimens  having  three  crack 
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orientations  with  respect  to  the  hole  and  many  different  crack  shapes. 
Two  hole  diameter  to  plate  thickness  ratios  were  duplicated  and  speci- 
mens with  load  transfer  to  a fastener  were  tested.  Agreement  between 
these  experiments  and  the  present  theory  at  specific  locations  on  the 
crack  periphery  was  good  in  most  cases.  This  provided  verification  of 
the  theory  for  crack  locations  with  respect  to  the  hole  and  plate 
loading  conditions  not  considered  in  previous  experiments. 

8.2  Reconwendations 

The  finite  element-alternating  method  presented  in  this  work  uses 
as  one  solution  that  of  stresses  in  an  infinite  body  due  to  pressure  on 
an  embedded  elliptical  crack  where  this  pressure  is  described  by  a ten- 
term  polynomial.  In  the  present  solution,  a least-square  fit  to  the 
actual  pressure  produces  the  required  polynomial  coefficients.  This 
polynomial  approximation  accurately  predicts  the  total  force  applied  to 
the  crack  and  correctly  reflects  trends  in  the  stress  gradients  yet 
subtleties  in  the  pressure  are  not  precisely  reproduced.  The  resulting 
magnification  factors  are  probably  slightly  in  error.  A more  accurate 
representation  of  this  portion  of  the  finite  element-alternating  method, 
such  as  a Green's  function  solution,  would  be  an  improvement. 
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Figure  1.2  Double  Through-Cracks  Emanating  From  a Hole  In  a Plate 
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(a)  Single  quarter-elliptical  corner  crack 


(b)  Double  quarter-el  1 Iptical  corner  cracks 
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(c)  single  scffll-ell Iptical  embedded  crack 


Figure  2.1  Crack  Types  and  Locations  Relative  to  a Fastener  Hole 
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Figure  2.6 


Megniflcation  Factor  Versus  Angle  from  Front  Surface 
Single  Corner  Cracks 
a/t-0.9,  d/t-1.0,  Q/P-0.0 


kuc-.i  r««r  rnos*  njffk'i 


Figure  2.7 


Magnification  Factor  Versus  Angle  from  Front  Surface 
Single  Corner  Cracks 
a/fO.75,  d/fl.O,  O/P-0.0 
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Figure  2.8  Magnification  Factor  Versus  Angle  from  Front  Surface 

Single  Corner  Cracks 
a/fO.5,  d/fl.O,  Q/P-0.0 
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figure  2,1 
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Figure  2.15  Co^rison  of  Magnification  Factors  at  the  Front  Surface,  Single  Corner  Crack,  Open  Hoi 
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Figure  3.4  Megnif Icatlon  Factor  Versus  Angle  front  Front  Surface 

Double  Corner  Cracks 
a/fO.9,  d/fl.O,  Q/P-0.0 
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figure  3.5 


M«gn1 f 1c8tion  Factor  Versus  Angle  from  Front  Surface 
Double  Corner  Cracks 
a/fO.75,  d/fl.O,  O/P-0.0 
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M«gn1f1c*t1on  Factor  Versus  Angle  from  Front  Surface 
Double  Corner  Cracks 
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PiR*«ET»IC  INCLE  rnow  fHONT  SUHFACt 


figure  3. 
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Magnification  Factor  Versus  Angle  from  Front  Surface 
Double  Corner  Cracks 
a/fO.SO,  a/t»0.?5.  «l/t*1.0,  Q/P-0.0 
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Fioure  3.9  Magnification  Factor  Versus  Angle  from  Front  Surface 

Double  Corner  Cracks 
a/c*0.25,  a/fO.125,  O/P-0.0 
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figure  ^-3 


Magnification  factor  Versus  Angle  from  Minor  Axis 
Embedded  Cracks 
a/f0.125.  d/fl.O,  0/P*0.0 
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Magnification  Factor  Versus  Angle  from  Minor  Axis 
Embedded  Cracks 
a/t*0.?5,  d/t«1.0,  0/P*0.0 
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Figure  4. 


Magi'ific4tion  Factor  V'^sus  Angle  from  Minor  Axis 
Embedded  Cracks 
*/c»0.5,  d/fl.O,  0/P*0.0 
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Magnification  Factor  Versus  Angle  from  Minor  Axis 
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Figure  5.2  Tangential  Stress  Along  the  Minimum  Section  In  a Plate 
Under  Uniaxial  Loading,  Open  Hole,  i*0.2 
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figure  5.7  »iagn1fic«t1on  Factor  Versus  Angle  froti'  Front  Surface 

Single  Corner  Cracks 
a/c*0.S,  a/t*0.25,  d/t*0.5 
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figure  5. 


Magnification  Factor  Versus  Angle  from  Front  Surface 
Single  Corner  Cracks 
a/c*2.5,  a/fO.85,  d/t*1.0 
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6.2  Experimental  Loading  Arrangements 
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I F1<)ure  6.3  Typtcal  Fractured  Surfaces  In  PMW 

* Specimens  with  Cracks  Emanating 
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Figure  6.6  PMMA  Baseline  Test,  Plate  No. 


Figure  6.7  Comparison  of  Theory  and  PMMA  Experiments 
Single  Corner  Crack,  d/t  = I.O,  Q/P  - 0.0 
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Fiqiirp  6.10  Compa'''son  of  Theory  and  PftW  Experiments 
Single  Corner  Crack,  d/t  = 1.0,  0/P  = 1.0 
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Double  Corner  Cracks,  d/t 
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Double  Corner  Cracks,  d/t 
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comparison  of  Theory  and  PftIA  Experiments 
Embedded  Crack,  d/t  = 1.0,  0/P  ^ 0.0 
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Figure  6.17  Contparison  of  Theory  and  PMMA  Experiments 
Single  Corner  Crack,  d/t  = O.S,  Q/P  = 0.0 


Fiqure  <i.lR  Comparison  of  Theory  and  PHMA  Experiments 
Sinqfe  Corner  Crack,  d/t  = 0.5,  0/P  = 0.0 


f<(jure  6.19  Comparison  of  Theory  and  f’WM  Experiments 
Double  Corner  Crack,  d/t  = 0.5,  Q/P  = 0.0 


160 


•*-c. 


UilD 

if 


o u 


161 


Figure  Comparison  of  Theory  and  PWW  Experiments 

Embedded  Crack,  d/t  = 0.5,  Q/P  = 0.0 
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Figure  6.23  Difference  Between  Theoretical  and  Experimental 
Results  as  a Function  of  Maximum  Local  Crack 
Plane  Normal  Stress  in  the  Experiments 


Parametric  Plot  of  Theory  to  Approximate  Slow  Growth 
Single  Corner  Crack,  c/t=0.S,  d/t=1.0 


Single  Corner  Crack,  d/t«1.0 


APPENDIX  A 


FINITE  ELEMENT-ALTERNATING  METHOD  CALCULATIONS 
AND  STRESS  FUNCTION  COEFFICIENTS 


Case 

Crack  Configuration 

d/t 

Q/P 

a/t 

a/c 

1 

Single  corner 

1 .0 

0.0 

0.5 

2.0 

2 

0.5 

1.5 

3 

0.5 

0.75 

4 

0.5 

0.5 

5 

0.5 

3.0 

6 

0.376 

0.446 

7 

0.9 

3.0 

8 

0.9 

2.0 

9 

0.9 

1.5 

10 

0.75 

2.0 

11 

0.75 

1 .5 

12 

0.25 

0.25 

13 

0.125 

0.25 

14 

0.2 

2.0 

15 

0.45 

0.9 

16 

0.55 

1.1 

17 

0.3 

1.4 

18 

0.4 

1.4 

19 

0.5 

1.4 

20 

0.6 

1.4 

21 

0.7 

1.4 
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Case 

Crack  Configuration 

d/t 

Q/P 

a/t 

a/c 

22 

Single  corner 

1.0 

0.0 

0.8 

1.4 

23 

0.85 

2.5 

24 

0.25 

0.5 

25 

0.1 

0.25 

0.5 

26 

0.5 

0.25 

0.5 

27 

1 .0 

0.25 

0.5 

28 

1 .0 

0.85 

2.5 

29 

1 .0 

0.2 

0.5 

30 

1 .0 

0.2 

2.0 

31 

1 .0 

0.3 

0.5 

32 

0.5 

0.0 

0.75 

3.0 

33 

0.5 

0.0 

0.25 

0.5 

34 

0.5 

1.0 

0.25 

0.5 

35 

Double  corner 

1.0 

0.0 

0.85 

2.5 

36 

0.125 

0.25 

37 

0.25 

0.5 

38 

0.5 

0.75 

39 

0.9 

2.0 

40 

0.75 

2.0 

41 

0.5 

2.0 

42 

0.5 

1 .5 

43 

0.75 

1.5 

44 

0.9 

1.5 

45 

0.15 

0.5 

46 

0.2 

0.5 

47 


0.1 


0.23 


0.6 
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Case 

Crack  Configuration 

d/t 

. Q/P 

a/t 

a/c 

48 

Double  corner 

1.0 

1 .0 

0.25 

0.6 

49 

1 .0 

0.85 

2.5 

50 

1 .0 

0.125 

0.25 

51 

1 .0 

0.25 

0.5 

52 

0.5 

0.0 

0.21 

0.75 

53 

0.5 

0.0 

0.75 

3.4 

54 

0.5 

1 .0 

0.25 

0.5 

55 

Embedded 

1 .0 

0.0 

0.125 

0.5 

56 

0.125 

0.25 

57 

0.125 

0.33 

58 

0.1875 

0.5 

59 

0.25 

0.67 

60 

0.25 

0.75 

61 

0.375 

0.75 

62 

0.25 

0.5 

63 

1.0 

0.125 

0.5 

64 

0.5 

0.0 

0.125 

0.5 

65 

0.5 

1.0 

0.125 

0.5 
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APPENDIX  B 


CRACK  OPENING  DISPLACEMENT  EQUATION 


I d b d b 


- r-'  ■ 

. C-o  f - 16xikv}^1  _ r r 

"'°LaV  aV  J L aV 

P r 24x(wv)V^  16x2(pv)^l\ 

L a^b**  ■ a«b^  Jj 


ley^CuvJ**' 


8y(pv)^'^^  16x^y(tiv)^ 


where:  wv  = a^b^  - b^x^  - a^y^ 

a = crack  semi-major  axis 
b * crack  semi -ml nor  axis 

X = crack  coordinate  in  major  axis  direction,  origin  at  crack 
center 

y = crack  coordinate  In  minor  axis  direction,  origin  at  crack 
center 

n * Polsslon's  ratio  (Stress  function  coefficients  In 
Appendix  A generated  with  n 0.25) 
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APPENDIX  C 


BENCHMARK  PROBLEM 

As  three-dimensional  fracture  mechanics  becomes  more  sophisticated 
and  capable  of  producing  usable  solutions,  the  need  arises  for  a 
standard  problem  configuration,  or  a so-called  "benchmark  problem," 
which  can  be  used  by  all  Investigators  In  the  field.  Such  a problem 
statement  Is  necessary  so  that  solutions  obtained  by  the  many  varied 
methods  can  be  compared,  one  against  another,  without  the  uncertainties 
of  interpolation  or  extrapolation  of  results. 

The  Workshop  on  Three-Dimensional  Fracture  Analysis  held  in 
Columbus,  Ohio,  April  26-28,  1976,  and  attended  by  most  of  the  leading 
three-dimensional  fracture  mechanics  investigators  * United  States, 
defined  several  of  these  "benchmark  problems."  For  ti...  case  of  a flaw 
adjacent  to  a fastener  hole  in  an  infinite  plate,  the  following  config- 
uration was  agreed  upon: 

Single  corner  flaw 
a/c  = 2.0 
a/t  = 0.2 
d/t  = 1.0 

The  method  described  earlier  in  this  paper  was  applied  to  the  above 
stated  "benchmark  problem"  and  results  are  shown  in  the  figures  that 
follow. 

For  the  open  fastener  hole  case,  with  remote  uniaxial  plate  load- 
ing, Figure  C-1  shows  plotted  values  of  magnification  factor  versus 
parametric  angle  from  the  front  surface  of  the  plate  while  Table  C-1 
lists  these  same  values  in  tabular  form  for  each  ten  degrees  of 
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parametric  angle.  Similarly,  Figure  C-2  and  Table  C-1  show  results  for 
the  case  of  all  the  remote  load  transferred  to  the  fastener. 

Several  interesting  observations  can  be  made  regarding  the 
"benchmark  problem"  solutions  presented  in  this  section.  The  open 
fastener  hole  case  displays  a higher  stress  intensity  factor  at  the 
front  surface  and  intermediate  locations  than  at  the  hole.  This  same 
trend  was  also  noted  earlier  in  this  work  with  other  flaws  having 
high  a/c  and  low  a/t  ratios.  Figure  C-3  is  a parametric  plot  of  Kj/a 
versus  a/t  for  a/c  = 2.0  and  reversal  in  the  location  of  maximum  stress 
intensity  is  clearly  seen  to  occur  at  an  a/t  ratio  very  close  to  that 
of  the  "benchmark  problem."  Of  perhaps  more  significance  is  a similar 
trend  seen  in  the  problem  of  all  the  remote  load  transferred  to  the 
fastener.  Since  the  crack  dimension  is  so  small  along  the  front  sur- 
face and  the  crack  shape  is  slim,  the  entire  crack  lies  in  an  area  of 
elevated  stress.  Thus,  the  crack  does  not  extend  out  into  an  area 
where  the  steep  stress  gradient  would  cause  a dramatic  reduction  in  the 
stress  intensity  factor  in  the  vicinity  of  the  front  surface.  With 
respect  to  the  "benchmark  problem"  solutions  presented  in  this  section, 
should  a/c  be  slightly  decreased  (holding  a/t  constant)  or  a/t  slightly 
increased  (holding  a/c  constant),  it  seems  apparent  that  the  location 
of  maximum  stress  intensity  would  shift  to  the  hole  intersection. 
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TABLE  C-1 

BENCHMARK  PROBLEM  RESULTS 


Open  Fastener  Hole 

100%  Load  Transfer 
to  Fastener  Hole 

Parametric  Angle  From 
Front  Surface 

a/iri’ 

0 

1 .2684 

0.4073 

10 

1.3223 

0.4211 

20 

1.3450 

0.4244 

30 

1.3468 

0.4210 

40 

1.3353 

0.4138 

50 

1.3133 

0.4043 

60 

1.2799 

0.3923 

70 

1 .2357 

0.3780 

80 

1.1897 

0.3638 

90 

1.1574 

0.3544 

s 
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